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Firms  are  increasingly  dependent  on  data  networks  to  complete  day-to-day 
operations  such  as  video  conferencing,  chat/online  customer  services,  real  time 
browsing  /E-Commerce  applications  and  other  data  applications.  In  general,  a firm 
performs  two  types  of  tasks  over  a network,  time-fixed  and  size-fixed.  These  tasks  might 
have  to  be  completed  using  networks  with  different  quality  of  service  levels  (QoS)  and 
pricing. 

We  investigate  the  customer’s  optimal  behavior  when  there  are  quality  related 
costs  in  addition  to  the  cost  of  acquisition  while  using  networks  with  bandwidth  and  QoS 
constraints.  We  consider  the  customer’s  cost  minimization  problem  assuming  that  an 
environment  exists  in  which  network  capacity  can  be  purchased  at  competitive  prices 
from  different  suppliers  with  different  service  quality.  We  model  QoS  as  a function  of  loss 
probability,  delay  and  jitter.  We  also  assume  that  the  suppliers  can  offer  any  quality  of 
service  and  some  capacity  at  competitive  prices.  Therefore,  the  customer  has  access  to 
sufficient  resources  with  different  QoS  levels,  capacities,  and  prices  to  perform  tasks 


IX 


with  different  QoS  and  capacity  requirements.  Existing  literature  on  QoS  mostly  focuses 
on  supplier  side  issues.  We  focus  on  the  customer’s  point  of  view  at  a tactical  level.  This 
study,  to  the  best  our  knowledge,  is  among  the  first  to  formulate  and  solve  this  problem. 

First,  we  analyze  the  problem  with  all-you-can-send  pricing  where  the  supplier 
charges  a fixed  price  for  a specific  bandwidth  and  duration.  We  show  that  the  resulting 
problem  is  at  least  as  hard  as  the  two-dimensional  bin-packing  problem.  We  develop  a 
heuristic  (Heuristic  A)  that  can  take  advantage  of  the  trade-off  between  the  cost  of 
acquiring  resources  and  the  opportunity  cost  of  degradation  in  realized  quality.  We 
analyze  how  a customer  behaves  as  a decision  maker  under  different  pricing  settings 
through  experiments.  The  customer’s  strategy  varies  based  on  prices,  capacity  and 
quality  of  resources.  Given  the  opportunity,  customers  mix  and  match  resources  to 
optimize  their  position. 

The  underlying  problem  is  a computationally  intractable  problem.  A relaxation  of 
this  problem  can  actually  create  workable  solutions  for  real  life  situations.  Qn  that  end 
we  develop  a lower  bound  (LB)  formulation  and  solve  it  optimally.  The  solutions  are 
reasonable  when  prices  set  randomly. 

A second  formulation,  which  is  a more  realistic  relaxation  of  the  original  problem, 
is  solved  optimally  by  the  Generalized  Bender’s  Decomposition  (GBD)  method.  We 
show  that  GBD  solves  the  small  problems  in  reasonable  time  and  the  solutions  are  very 
realistic.  However,  more  work  is  needed  to  make  it  work  in  large  problems. 

Finally  we  analyze  a different  formulation  using  a pricing  scheme  called  tax-band 
pricing  in  which  the  total  price  of  a resource  is  a convex  piecewise  linear  function  of 
capacity.  We  show  that  this  problem  is  easier  to  solve. 
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CHAPTER  1 
INTRODUCTION 

Firms  today  use  telecommunications  in  a variety  of  ways.  In  particular,  they  are 
increasingly  dependent  on  data  networks  to  complete  routine  day-to-day  operations. 
Carrying  out  these  daily  activities  requires  completing  a set  of  tasks  whose  demands  on 
the  network  are  quite  different  creating  a heterogeneous  traffic  stream  on  the  underlying 
network.  Typical  tasks  require  running  real  time  streaming  applications  such  as  video 
conferencing,  voice  over  TCP/IP,  real  time  transactional  applications  such  as  online 
procurement  applications,  asynchronous  messaging  applications  such  as  email,  and 
other  batch  transactional  applications.  In  order  to  carry  out  these  applications,  firms  rely 
on  data  networks.  For  example,  traditional  face-to-face  conferencing  is  preferred  less, 
while  virtual  conferencing  is  becoming  popular  increasing  the  reliance  on  data  networks 
(Figure  1-1). 
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Figure  1-1.  The  increase  in  teleconferencing  usage 
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In  a survey  carried  out  by  MCI,  one  of  the  American  telecommunications  firms, 
69%  of  the  people  preferred  a virtual  meeting  to  traveling  because  it  saved  time,  while 
only  37%  preferred  a virtual  meeting  since  it  saves  money  (Being,  2004).  By  using 
teleconferencing  services  one  can  convene  virtual  meetings  attended  online  by  all  the 
right  people,  wherever  they  are  and  convey  strategic  information  quickly,  easily,  and 
cost-effectively  (AT&T  Alascom,  2004).  Moreover,  it  reduces  travel  time  and  costs,  and 
increases  the  productivity  of  workers  who  participate  fully  without  leaving  their  desks. 

Each  application  that  relies  on  data  networks  has  different  capacity  and  quality  of 
service  (QoS)  requirements  with  different  sensitivity  to  due-date  constraints.  Each  is 
affected  differently  by  the  variability  in  the  network  service  availability  and  speed.  This 
poses  challenging  strategic  and  tactical  problems  on  the  firm  as  a user  of  a network  in 
an  environment  where  network  resources  are  leased  from  third  parties.  Even  though  the 
real  time  mechanics  are  quite  involved  and  hard  to  determine  in  advance  we  believe,  for 
a sound  telecommunications  policy,  the  firm  (customer)  has  to  determine  which 
resources  to  acquire  and  how  to  assign  existing  tasks  to  these  resources.  Such  an  effort 
requires  solving  a cost  minimization  (telecommunications  and  opportunity  costs)  problem 
subject  to  a variety  of  capacity,  due-date  and  quality  constraints  based  on  known  tasks 
to  be  completed  during  a single  planning  period. 

A network  operator  (supplier)  faces  a heterogeneous  set  of  tasks  with  different 
demands  on  capacity  and  quality.  The  network  operator  has  to  manage  its  resources 
and  pricing  such  that  it  can  maximize  its  profits  subject  to  capacity  requirements  and 
potential  competition  from  other  network  operators.  Today’s  networks  are  designed  on 
the  assumption  that  users  behave  cooperatively;  in  other  words  they  reduce  their 
transmission  rate  when  they  detect  congestion.  If  the  users  do  not  reduce  their 
transmission  rate  fast  enough,  the  congestion  causes  collapse  of  the  network.  It  is 
important  to  control  the  network  traffic  instead  of  relying  on  users  and  waiting  for  them  to 
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reduce  their  transmission  rate.  One  solution  for  control  is  charging  the  end  users  for  the 
network  use  (Rasmusson,  2002).  Network  operators  can  control  network  congestion  and 
divide  the  network  resources  among  users  by  using  charges.  Charging  generates 
income,  feedback  and  control  so  that  network  performance  could  be  improved 
(Courcoubetis  et  al.  2000).  Given  services,  prices  and  charging  methods,  the  operator  of 
the  network  determines  the  quantity  of  services  (segmentation)  knowing  that  the  users’ 
demand  fluctuates.  Network  providers  charge  higher  prices  for  real  time  services  and 
low  prices  for  the  best  effort  services  (Courcoubetis  et  al.  2000). 

Suppliers  can  offer  different  charging  schemes  for  the  capacity  they  sell.  For 
example,  ISPs  offer  a combination  of  fixed  price  for  a fixed  maximum  bandwidth  (i.e.. 
Cable,  DSL)  or  pay  per  hour  (or  minute)  for  again  a maximum  bandwidth.  Wireless 
phone  companies  charge  for  text  messaging  based  on  bytes  sent.  Certain  calling  plans 
offer  different  per  minute  charges  for  phone  conversations  depending  on  when  the  call  is 
placed. 

The  strategy  of  the  customer  will  vary  based  on  the  prices  and  the  associated 
quality  offered  by  the  suppliers  in  the  marketplace.  Given  the  opportunity,  the  customer 
will  mix  and  match  resources  to  optimize  its  position.  For  example,  while  resources 
purchased  for  a fixed  price  can  be  used  for  data  applications  with  lower  QoS 
requirements,  more  expensive  ones  might  be  utilized  for  real  time  applications  with  high 
QoS  requirements.  In  doing  so,  the  firm  can  acquire  resources  with  enough  QoS  levels 
and  incur  less  costs  due  to  degradation  in  the  quality  of  certain  tasks,  which  we  call 
opportunity  costs.  The  firm  can  acquire  fewer  resources  and  incur  less 
telecommunication  cost  by  decreasing  QoS  level  of  certain  tasks  and  incurring  higher 
opportunity  cost.  Therefore,  the  problem  that  the  firm  has  to  solve  is  a cost  minimization 
problem  that  reflects  a trade  off  between  cost  of  acquiring  resources  and  the  opportunity 
cost  of  degradation  in  realized  quality. 
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It  is  in  this  setting  we  wish  to  nnodel  the  problem  from  the  customer’s  point  of 
view  at  a tactical  level.  We  consider  the  customer’s  cost  minimization  problem  assuming 
that  an  environment  exists  in  which  the  firm  can  buy  network  capacity  at  competitive 
prices  from  different  suppliers  with  different  service  quality.  We  also  assume  that  the 
suppliers  can  offer  any  quality  of  service  and  some  capacity  at  competitive  prices. 

Bandwidth  Exchange  is  a network  bandwidth  market  in  which  companies  buy 
and  sell  capacity.  It  is  managed  by  broker  companies  such  as  TFS  Telecom,  Band-X, 
Amerex  Bandwidth,  Arbinet,  and  RateXchange  (Rasmusson,  2002).  According  to 
RateXchange,  the  bandwidth  trading  market  is  expected  to  grow  to  $12  billions  within 
five  years  (Telecommunications  Industry  Association,  2000).  Bandwidth  broker  is  a 
neutral  trading  floor  for  providers  to  buy  and  sell  capacity  on  a month-to-month  basis, 
called  bandwidth  on  demand  (Wagner,  2001).  It  basically  acts  as  a broker  between 
network  carriers  (backbone  providers)  such  as  AT&T,  Global  Crossing,  Qwest 
Communications  International,  Williams  Communications  (Ames,  2001)  and  the 
worldwide  bandwidth  buying  community.  In  order  to  complete  the  tasks  required  for  daily 
activities  the  customer  has  to  first  acquire  (basically  lease)  the  capacity  (bandwidth)  from 
either  broker  firms  or  backbone  providers  (i.e.,  MCI,  Sprint,  UUNET  etc.),  then  allocate 
the  tasks  into  that  capacity  (Figure  1-2). 

The  quality  of  service  (or  the  lack  of  it)  provided  by  the  supplier  adversely  affects 
the  customer  in  two  ways.  First,  certain  tasks  might  be  delayed  beyond  planned 
deadlines  creating  unacceptable  schedules.  Second,  the  realized  quality  of  a task  such 
as  a videoconference  might  be  less  than  acceptable  creating  an  opportunity  cost  for  the 
customer.  Each  task  has  different  importance  and  therefore  a different  opportunity  cost. 
Customers  want  to  designate  different  types  of  traffic  for  different  types  of  treatments 
(Carpenter  and  Nichols,  2002). 
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In  this  study,  we  consider  the  customer  as  a rational  decision  maker  who  has  to 
make  two  decisions;  (1)  to  determine  which  resources  to  acquire,  and  (2)  how  to  assign 
existing  tasks  to  these  resources.  Existing  research  focuses  only  on  infrastructure  and 
pricing  issues.  There  is  no  treatment  of  this  problem  as  a resource  selection  (vendor 
selection)  problem  in  the  literature.  From  a managerial  point  of  view,  solving  this  problem 
will  help  us  to  understand  how  a customer  behaves  under  different  pricing  conditions 
with  different  QoS  and  capacity  constraints.  Moreover,  we  wish  to  understand  the 
properties  of  an  optimal  resource  allocation  strategy.  The  main  theoretical  objective  of 
this  study  is  to  develop  efficient  optimal  algorithms  or  heuristics. 

The  rest  of  this  dissertation  document  is  organized  as  follows.  In  Chapter  2,  a 
review  of  quality  of  service  issue  for  telecommunication  networks  is  given.  A literature 
review  of  quality  of  service  and  pricing  issue  for  telecommunication  networks  is  given  in 
Chapter  3.  The  formulations  of  the  main  problems  are  given  in  Chapter  4.  Chapter  5 
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focuses  on  a specific  formulation  that  uses  all  you  can  send  charging  scheme.  A lower 
bound  procedure  using  the  Generalized  Benders  Decomposition  method  is  given  in 
Chapter  6.  A cost  minimization  problem  with  the  tax-band  charging  scheme  is  given  in 
Chapter  7.  Summary  and  future  research  directions  are  presented  in  Chapter  8. 


CHAPTER  2 

QUALITY  OF  SERVICE  FOR  TELECOMMUNICATION  NETWORKS  OVERVIEW 

Quality  of  service  is  a very  popular,  overloaded  term  and  is  looked  at  from 
different  perspectives  by  the  networking  and  the  application-development  communities 
(Miras,  2002).  In  networking,  “Quality  of  Service”  refers  to  the  ability  to  provide  different 
treatment  to  different  classes  of  traffic.  The  primary  goal  is  to  increase  the  overall  utility 
of  the  network  by  granting  priority  to  higher-value  or  more  performance-sensitive  traffic. 
“Priority”  means  either  reduced  loss  probability  or  preferential  queuing  at  congested 
interfaces  (Miras,  2002). 

The  view  of  QoS  that  application  developers  and  application  users  have  is  more 
subjective.  It  is  seen  as  something  that  will  improve  application  performance.  However, 
QoS  may  or  may  not  improve  an  individual  application’s  performance.  Results  are  highly 
dependent  on  the  relationship  between  a particular  application’s  realized  quality  (that  is, 
how  pleasant  or  unpleasant  the  experience  is  to  the  user  of  the  application,  e.g.,  visual 
quality  of  a displayed  video  sequence),  and  the  network  performance  it  experiences. 
Realized  quality  may  reflect  the  application’s  ability  to  perform  its  task  (for  example,  in  IP 
telephony,  whether  or  not  good  conversation  is  achieved)  or  generate  user  interest.  The 
application’s  performance  is  also  dependent  on  whether  a particular  individual  can  afford 
the  extra  cost  of  priority  treatment  (Miras,  2002).  In  many  Internet  service  markets,  the 
cost  of  such  priority  treatment  exceeds  the  cost  of  upgrading  to  faster  best-effort  service. 

The  QoS  measures  that  apply  to  the  network  do  not  necessarily  coincide  with  the 
user’s  perception  of  QoS  parameters.  A user  concerned  about  an  application’s 
performance  is  not  necessarily  concerned  with  the  details  of  how  a network  service  is 
implemented  and  how  it  performs.  Application  performance  deterioration  should  be 
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expressed  in  terms  that  focus  on  user-perceivable  effects  rather  than  on  their  origin 
within  the  application  end-to-end  path  (Miras,  2002).  For  example,  the  end-user  of  an 
application  sees  only  latency,  and  cannot  distinguish  the  cause  of  it,  whether  it  comes 
from  the  network  or  from  the  processing  done  on  the  end-system.  Thus,  from  a user 
perspective,  network  QoS  performance  is  hidden  in  the  end-to-end  application-level 
performance.  Application-level  performance  measures  include  throughput,  latency  and 
service  availability. 

Several  other  factors  in  the  end-to-end  application  path  may  result  in  undesired 
changes  of  performance  parameters:  the  operating  system’s  inability  to  support  the 
application,  continuity  of  service,  erroneous  application  and  protocol-stack 
implementations,  the  usage  environment  (e.g.,  faulty  equipment),  etc.  Given  that  such 
an  interaction  exists  among  the  application  software,  the  hardware  and  the  network  the 
only  objective  view  of  QoS  parameters  are  those  that  are  controlled  by  the  network 
provider. 

2.1  QoS  Performance  Measures 

Throughput  and  latency  are  more  important  factors  for  the  user  point  of  view 
performance  measures.  Availability  and  continuity  of  service  describe  the  requirement 
for  uninterrupted  service  with  acceptable  quality.  Network  availability  is  a reliability  issue 
and  we  will  assume  that  service  is  available  at  all  times. 

2,1.1  Throughput 

Bandwidth  defines  the  available  data  speed  within  the  network.  It  represents  the 
portion  of  the  available  capacity  of  an  end-to-end  network  path  that  is  accessible  to  the 
application  or  data  flow.  Consequently,  the  numbers  of  bits  that  are  injected  into  the 
network  by  the  various  flows  of  an  application  have  to  be  adjusted  accordingly. 
Bandwidth  is  limited  by  physical  media  connection  rates  and  network  congestion. 
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Throughput  is  a measure  of  the  rate  at  which  data  can  be  sent  through  the 
network,  representing  the  available  network  capacity  in  bits  per  second.  Throughput  is 
the  effective  share  of  bandwidth  that  the  application  receives  from  the  network. 
Bandwidth  and  throughput  are  used  interchangeably.  In  fact,  bandwidth  is  the 
throughput  capability  of  the  underlying  hardware,  which  poses  an  upper  bound  on 
throughput  (Comer,  2001). 

Throughput  affects  the  realized  quality  and  effective  transmission  rate  so  that  it 
impacts  latency.  Non-real  time  traffic  guarantees  minimum  throughput  and  utilizes  all 
available  throughput  at  any  time  (best  effort).  However,  real  time  traffic  requires 
guaranteed  constant  data  rate  and  delay. 

2.1.2  Latency 

Latency  is  the  end-to-end  delay  that  an  application  experiences.  It  has  a direct 
impact  on  user  satisfaction.  Latency  can  be  defined  as  a function  of  packet  loss 
probability,  delay  and  jitter. 

2. 1.2.1  Packet  loss 

Packet  loss  represents  data  dropped  or  irrecoverably  damaged  by  the  network. 
Packet  loss  is  typically  the  result  of  data  collision  and  buffer  overflows,  usually  due  to 
excessive  congestion  in  the  network.  Packet  loss  is  defined  as  the  fraction  (or 
percentage)  of  IP  data  packets  lost  somewhere  along  the  path  from  the  source  to  the 
destination  (Miras,  2002).  Packet  loss  is  a fundamental  problem  in  computer  networks 
because  transmission  errors  can  corrupt  bits,  making  the  frame  invalid.  When  a receiver 
detects  such  a problem,  it  discards  the  frame.  Depending  on  the  type  of  application,  the 
receiver  can  ignore  the  lost  frames  or  send  a retransmission  request  to  the  sender.  Most 
data  applications  are  extremely  sensitive  to  lost  packets.  A financial  statement,  report,  or 
electronic  funds  transfer  is  seriously  affected  if  data  are  missing  or  corrupted.  Data 
networks  must  have  an  ability  to  detect,  correct,  or  retransmit  corrupted  or  lost 
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information.  It  is  very  important  to  detect,  correct,  or  retransmit  corrupted  or  lost 
information  within  data  networks  (Ragsdale  et  al.  2000).  For  Voice-over-IP  (VoIP) 
networks,  voice  packets  have  a high  tolerance  for  truly  random  packet  loss.  VoIP 
applications  can  reconstruct  the  packets  if  the  number  of  lost  packets  is  minimal.  As  a 
rule  of  thumb,  the  packet  loss  rate  should  not  be  more  than  5%  in  VoIP  networks. 
Otherwise,  the  voice  quality  will  be  compromised  and  there  will  be  a significant  loss  in 
voice  quality  (Ragsdale  et  al.  2000;  Transition  Networks,  2003).  The  loss  pattern  or  the 
loss  period  is  also  equally  important.  The  most  notable  observation  concerning  loss 
patterns  in  today’s  Internet  is  that  losses  appear  in  bursts  that  are  usually  short,  but 
occasionally  lengthier  when  there  is  extreme  congestion.  If  loss  ratio  exceeds  the 
service  level  agreement  (SLA)  loss  ratio  for  the  user,  a possible  SLA  violation  is  reported 
(Habib  et  al.  2003).  Packet  loss  can  be  overcome  with  silence  substitution,  noise 
substitution,  packet  repetition,  packet  interpolation,  and  frame  interleaving. 

For  real  time  applications,  the  user-perception  of  information  loss  does  not 
necessarily  coincide  with  data  packet  loss  at  the  network  level.  It  might  be  data  loss  at 
the  application  level  (for  example,  loss  of  signal  fidelity  due  to  encoding),  or  it  might  be 
completely  user-specific.  It  is  described  as  a “user  perceivable”  loss  of  information  that 
leads  to  user  discomfort,  lack  of  timely  judgment  and  reaction  needed  to  execute  the 
task  successfully,  annoyance,  disorientation,  lack  of  interest,  etc.  Therefore,  information 
loss  corresponds  to  the  amount  of  information  (visual  content  as  the  user  sees  it,  loss  of 
audio  clarity,  etc.)  that  a user  perceives  as  missing. 

2.1. 2.2  Delay 

Delay  specifies  how  long  it  takes  for  a data  unit  to  travel  across  the  network  from 
one  computer  to  another  (from  source  to  destination)  measured  in  seconds  or  fractions 
of  seconds  (Comer,  2001 ).  Network  delay  is  caused  by  a combination  of  network 
propagation  delay,  processing  delays  and  variable  queuing  delays  at  the  intermediate 
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routers  on  the  path  to  the  destination  host  (Miras,  2002).  Queuing  introduces  the  most 
delay  into  a data  network.  Each  network  element  performs  some  degree  of  queuing  to 
allow  the  device’s  processor  to  make  routing  or  switching  decisions.  Queuing  delays 
increase  during  periods  of  congestion,  since  the  processors  in  the  network  elements 
must  work  harder  to  handle  the  increased  volume  (Ragsdale  et  al.  2000).  Video 
transmission  is  less  sensitive  to  delay  and  data  loss.  However,  it  requires  more 
bandwidth,  and  synchronization  with  audio  is  an  important  issue.  Delay  is  an  important 
factor  for  quality  of  audio  transmission.  In  Table  2-1 , the  limits  of  delay  for  the  quality  are 
given  (Teitelbaum  and  Sadagic,  2002). 

Table  2-1 . Delay  guidelines  for  VoIP 


Qne-way  Delay 

Effect  on  Realized  Quality 

<100-150  ms 

Delay  not  detectable 

1 50-200  ms 

Acceptable  quality;  slight  delay  or  hesitation  noticeable 

Qver  200-300  ms 

Unacceptable  delay:  normal  conversation  impossible 

Delay  is  closely  tied  to  bandwidth.  Bandwidth  impacts  delay  by  limiting  the 
maximum  transfer  unit  (MTU)  size  for  a given  link.  Lower  speed  links  typically  have 
smaller  MTU  sizes,  since  the  cost  of  retransmitting  a large  packet  is  significant  on  a low 
speed  link.  Smaller  MTU  sizes  translate  to  more  packets  per  payload  byte,  which  in  turn 
translates  into  more  overhead  consuming  the  link’s  bandwidth  and  increased  delay.  Due 
to  more  packets,  more  overheads  consume  the  link’s  bandwidth  and  increases  delay 
(Ragsdale  et  al.  2000). 

The  primary  QoS  challenges  for  Voice-over-IP  are  delay  and  jitter.  Although 
individuals  may  disagree  on  how  much  delay  is  tolerable  in  an  audio  conversation,  it  is 
generally  accepted  that  150  ms  is  the  break  point  where  the  delay  becomes  annoying  for 
most  people.  Delay  is  not  that  critical  in  data  networks.  Waiting  for  a web  page  to  load  or 
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a transaction  to  complete  is  not  as  irritating  as  a silence  when  you  are  in  the  middle  of 
an  important  conversation.  While  it  is  acceptable  to  wait  5 seconds  loading  a web  page, 
it  is  impossible  to  tolerate  such  a delay  during  a conference  call  with  the  customer. 

2.1. 2.3  Jitter 

Jitter  represents  the  variance  in  delay.  Congestion  adds  to  the  delay  of  a network 
by  holding  packets  in  queues.  Variations  in  traffic  intensity  cause  variations  in 
congestion  and  in  queuing  delay  creating  jitter.  Delay  variation  is  usually  caused  by  the 
buffers  built  up  on  routers  during  periods  of  increased  traffic,  and  less  often  by  changes 
of  routing  due  to  failures  or  routing  table  updates.  The  audio  and  video  applications  are 
quite  sensitive  to  delay  and  jitter  while  data  services  are  insensitive  to  both  of  them 
(Ragsdale  et  al.  2000). 

Voice  traffic  cannot  tolerate  much  jitter.  Without  compensation,  jitter  has  the 
effect  of  speeding  up  or  slowing-down  audio  and  quality  of  the  conversation.  Jitter 
appears  to  the  user  as  awkward  pauses  in  the  middle  of  a word.  In  the  extreme,  jitter 
may  either  overrun  or  starve  the  audio  decoder,  leading  to  a loss  of  synchronization  with 
further  disruption  to  an  otherwise  orderly  conversation  (Transition  Networks,  2003).  The 
limits  of  jitter  for  the  quality  are  given  in  Table  2-2  (Miras,  2002). 

Table  2-2.  Jitter  guidelines  for  VoIP 


Jitter 

Effect  on  Realized  Quality 

<40  ms 

Jitter  not  detectable 

40-75  ms 

Good  quality,  but  occasional  delay  or  jumble  noticeable 

Over  75  ms 

Unacceptable:  too  much  jumble  in  the  conversation 

A jitter  buffer  added  to  the  Voice-over-IP  connection  is  a common  solution  to 
jitter.  Jitter  buffers  store  voice  packets  in  memory  for  a time  equal  to  the  average  delay 
variation.  Large  delay  variations  equate  to  large  jitter  buffers  and  large  built-in  delays. 
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Jitter  buffers  mitigate  delay  variation,  but  they  compound  delay  problems.  Consequently, 
it  is  necessary  to  keep  delay  jitter  to  a minimum  for  reasons  of  quality  and  cost. 

Although  certain  levels  of  jitter  can  be  absorbed  by  the  application  through  the 
use  of  buffering,  application  data  packets  may  reach  the  receiver  too  late  to  be  useful. 
Delay  variation  can  also  produce  temporal  inconsistency  in  presenting  the  data  to  the 
user  (e.g.,  video  freezing  and  then  starting  again)  that  results  in  the  user  becoming 
irritated  and  losing  the  sense  of  presence.  In  time-critical  operations  as  well  as  human- 
to-human  interactions,  this  is  worse  than  having  a constant  delay. 

Prioritizing  traffic  can  eliminate  delay  jitter.  However,  processing  must  be  fast  and 
packet  queuing  kept  to  a minimum  if  low  delay  and  jitter  are  to  be  achieved.  Othen/vise, 
the  prioritization  process  adds  to  the  delay  rather  than  mitigating  it  (Ragsdale  et  al. 

2000).  Jitter  and  loss  probability  are  related.  Sudden  changes  in  jitter  may  cause  loss 
(Teitelbaum  and  Sadagic,  2002). 

In  summary,  bandwidth  (a  proxy  for  throughput)  and  latency  will  be  modeled  in 
this  study.  Substituting  bandwidth  for  throughput  allows  us  to  decouple  latency  and 
throughput.  Raw  bandwidth  is  the  assumed  capacity  per  second  that  we  can  model 
explicitly.  Since  delay  and  jitter  require  packet  level  treatment  and  loss  probability  can  be 
incorporated  at  the  task  level,  we  will  model  loss  probability  explicitly  and  delay  and  jitter 
as  simple  levels  that  must  be  met.  We  show  that  packet  loss  reduces  effective 
bandwidth  and  causes  either  delay  or  quality  degradation. 

2.2  Traffic  Classification  and  Prioritization 

There  is  a widely  held  belief  that  advanced  applications  cannot  be  entirely 
accommodated  by  today’s  Internet,  and  that  it  is  necessary  to  have  a service  model  that 
offers  QoS  guarantees  to  flows  that  need  them.  There  is  another  belief  that  QoS  needs 
of  applications  can  be  sufficiently  met  by  an  over-provisioned  best  effort  network, 
combined  with  application  intelligence  to  adapt  to  the  changing  availability  of  network 
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resources  and  to  tolerate  loss  and  jitter.  Both  of  the  approaches  have  merits  and 
disadvantages.  It  is  probably  true  that  an  efficient  solution  lies  somewhere  between 
these  two  positions  and  favors  some  form  of  traffic  differentiation  (Miras,  2002).  It  is 
apparent  that  without  any  form  of  traffic  classification  and  prioritization,  network 
congestion  will  become  a problem,  affecting  QoS-sensitive  flows  and  reducing  the 
quality  of  the  corresponding  applications.  However,  the  selection  of  a suitable  network 
model  is  a complicated  function  of  several  factors,  such  as  the  criticality  of  the 
applications,  the  complexity  and  scalability  of  the  solution,  and  the  economic  model  or 
the  market  needs  (Miras,  2002). 

2.2.1  Prioritization 

Bandwidth  is  more  available  and  less  expensive  than  ever  before  due  to  the 
many  innovations  in  fiber  optics,  dense  wave  division  multiplexing  (DWDM),  and  high- 
speed electronics.  Even  so,  infinite  bandwidth  equates  to  infinite  cost,  neither  of  which  is 
practical  within  commercial  networks  (Ragsdale  et  al.  2000).  Prioritization  of  traffic  is  a 
compromise  between  infinite  bandwidth,  packet  processing,  and  cost. 

The  fundamental  idea  behind  QoS  is  that  the  network  handles  traffic  by 
considering  the  type  of  the  application  and  data  being  sent  (Rasmusson,  2002).  For 
example,  a real-time  data  stream  should  be  transmitted  with  minimal  delay  and  loss, 
while  batch  file  transfer  may  get  sufficient  service  when  the  network  is  busy 
(Rasmusson,  2002).  Ideally,  unlimited  bandwidth  would  eliminate  jitter  by  ensuring  that 
all  network  connections,  routers,  and  switches  have  sufficient  bandwidth  to  support  the 
most  intense  traffic  loads.  However,  networks  will  get  congested  at  peak  times  even 
when  there  is  seemingly  ample  bandwidth  available.  When  the  packet  network  becomes 
congested,  data  packets  have  to  be  delayed  or  dropped.  Prioritization  determines  which 
data  are  delayed  by  using  service  classes. 
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Prioritization  of  IP  packets  is  a computationally  intensive  task.  Processing  must 
be  fast  and  packet  queuing  kept  to  a minimum  if  low  delay  and  jitter  are  to  be  achieved. 
Othen/vise,  the  prioritization  process  adds  to  the  delay  rather  than  mitigating  it. 

2.2.2  Service  Classes 

Service  classes  are  used  to  handle  the  multi-service  traffic  over  the  network. 
Internet  service  models  have  been  developed  by  using  service  classes.  Each  class 
corresponds  to  one  or  more  type  of  application  and  each  class  has  its  own  delay,  jitter, 
and  transmission  rate  parameters.  Table  2-3  gives  an  example  of  service  classes  that  is 
described  by  Zander  (2002). 

As  we  mentioned  earlier  each  service  has  its  own  QoS  requirements  and  each 
network  technology^  has  its  own  challenges  to  meet  these  requirements.  Different 
technologies  require  different  solutions  creating  a different  level  of  management 
complexity.  For  example,  a common  problem  with  bandwidth  utilization  in  WAN  links  is 
the  lack  of  efficient  mechanisms  for  dealing  with  packet  loss  on  connection  overloads. 
As  a result,  all  types  of  traffic  flowing  through  the  link  suffer  substantial  deterioration  of 
quality.  Robust  QoS  policies  have  been  developed  to  enhance  the  IP  backbone's 
efficiency  and  reliability  and  to  reduce  packet  drops  and  retransmission  (Agilent 
Technologies,  2002).  The  Cisco  ITP  gateway,  for  example,  will  support  IP  Precedence 
and  differential  services  with  eight  QoS  classes  with  preferential  treatment  over  the  IP 
backbone,  thus  ensuring  predictable  SS7^  (Common  Channel  Signaling  System  No.  7) 
service  delivery  for  various  types  of  SS7  traffic  (Agilent  Technologies,  2002).  Real  time 
or  critical  applications  requiring  stable  network  bandwidth  can  be  assigned  higher 

' Here  we  use  technology  to  mean  the  underlying  medium,  protocol  and  topology  that  make  up 
the  network. 

2 

SS7  is  a standard  for  telecommunications  defined  by  the  International  Telecommunications 
Union  (ITU)  Telecommunication  Standardization  Sector  (ITU-T).  It  defines  procedures  and 
protocols  where  network  elements  in  the  public  switched  telephone  network  (PSTN)  exchange 
information  over  a digital  signaling  network  to  effect  call  setup,  routing  and  control. 
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priority,  which  guarantees  that  in  the  event  of  network  congestion  their  packets  will  not 
be  dropped  (lost),  and  assigned  bandwidth  parameters  and  QoS  will  be  preserved. 
Table  2-3.  Service  classes 


Service  Class 

Typical  applications 

Service  functional  characteristics 

Conversational 
Real  Time  (RT) 

Voice 

Preserve  time  relations  between  entities 
Stringent  preservation  of  conversational 
patterns  (low  delay) 

Streaming  RT 

Video/Audio 

streams 

Preserve  time  relations  between  entities 

Interactive 
Best  effort  (BE) 

Web-browsing 

Request-response  pattern 
Preserve  payload  (low  error  rate) 

Background 

BE 

File  transfer.  E-mail 

Not  time  critical 

Preserve  payload  (low  error  rate) 

A main  issue  in  the  QoS  is  the  bandwidth  management  that  directly  affects  the 
capacity  allocated  to  a task.  In  practice,  to  guarantee  the  flow  of  traffic,  allowable 
transmission  rates  have  been  defined.  Committed  information  rate  (CIR)  is  the  minimum 
amount  of  bandwidth  guaranteed  for  an  Ethernet  service.  Peak  information  rate  (PIR)  is 
the  maximum  rate  at  which  packets  are  allowed  to  burst  over  the  CIR  (Khandekar, 
2002).  Maximum  buffer  size  (MBS)  defines  the  maximum  offered  load  without  discard. 
MBS  specifies  the  maximum  amount  of  buffers  allocated  for  service  traffic  over  PIR.  If 
transmission  rate  is  greater  than  PIR  and  MBS  is  exceeded,  transmission  could  not  be 
processed.  If  transmission  rate  is  in  between  CIR  and  PIR  or  greater  than  PIR  but  MBS 
is  not  exceeded,  transmission  could  be  performed  but  delay  and  loss  probability  could 


increase. 
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2.3  Internet  Service  Models  for  Multi-Service  Traffic 

The  original  best  effort  LAN  protocols  were  designed  for  applications  such  as 
basic  connectivity  between  stations,  file  transfer,  e-mail  and  web  browsing.  The  best 
effort  strategy  is  to  simply  handle  packets  in  the  order  they  arrive,  and  discard  excess 
packets  that  do  not  fit  in  the  buffer  (Rasmusson,  2002).  In  the  best  effort  service,  the 
network  is  unable  to  guarantee  a delivery  deadline  time,  in  other  words  real-time 
guarantee,  because  it  is  impossible  to  know  how  much  retransmissions  is  necessary  due 
to  congestion  before  all  packets  arrive  the  destination  (Rasmusson,  2002). 

Multi-service  traffic  is  difficult  to  handle  efficiently  because  each  type  of  traffic 
requires  different  transfer  rate  and  each  has  a different  tolerance  for  delay  or  packet 
sequencing.  The  existing  Internet  service  (i.e.,  the  best  effort  service  of  IP)  cannot 
satisfy  the  QoS  requirements  of  emerging  multimedia  applications  (Li  et  al.  2000  and 
Transition  Networks,  2003).  To  support  multi-service  traffic  some  mechanisms  and 
service  models  have  been  proposed.  Integrated  service,  and  differentiated  service 
models  and  ATM  are  used  most  and  play  a primary  role  in  achieving  quality 
requirements. 

2.3.1  Integrated  Service  Model  (INT-SERV) 

Integrated  Services  (INT-SERV)  provides  a method  to  reduce  variability  in  delay 
and  reserving  a required  level  of  bandwidth  in  a network  to  support  QoS  effort  and  real- 
time protocols.  The  purpose  of  INT-SERV  is  to  guarantee  a service  level  for  delay 
sensitive  (real  time)  applications,  and  manage  loads  on  network  components  to  provide 
“enhanced  best-effort”  performance  (Li  et  al.  2000).  INT-SERV  has  four  components: 
flow  specification,  the  signaling  protocol,  admission  control,  and  packet  management 
(classification  and  scheduling).  INT-SERV  signaling  relies  on  Resource  Reservation 
Protocol  (RSVP)  to  make  service  class  reservations.  Two  classes  of  service  in  addition 
to  the  standard  “best  effort”  service  offered  by  packet  networks  are  defined  (Li  et  al. 
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2000  and  Transition  Networks,  2003).  The  guaranteed-service  (rate-sensitive)  class 
attempts  to  create  a virtual  circuit.  The  controlled-load  (delay-sensitive)  service  class 
attempts  to  control  the  delay  experience  by  packets.  For  a guaranteed-service  class,  a 
maximum  allowed  packet  loss  is  defined.  For  a controlled-load  class,  a maximum 
allowed  packet  delay  is  defined.  Both  classes  of  service  specify  the  parameters  of 
sustained  rate,  smallest  packet  size,  and  largest  packet  size.  The  fundamental 
assumption  of  INT-SERV  model  is  that  resources  (e.g.  bandwidth  and  buffer)  must  be 
explicitly  managed  for  each  real  time  applications  (Li  et  al.,  2000). 

RSVP  packets  travel  through  the  network  from  source  to  destination.  At  each 
router  along  the  path,  RSVP  requests  a service  class.  When  the  destination  receives  the 
packet,  it  returns  the  packet  to  the  source.  This  allows  the  source  to  know  the  minimum 
service  parameters  accepted  by  all  routers  in  the  path.  Routers  use  admission  control  to 
determine  whether  to  fully,  partially,  or  not  reserve  the  requested  bandwidth.  Admission 
control  is  largely  based  on  existing  commitments  and  available  bandwidth.  Routers 
classify  received  packets  to  determine  if  the  packet  matches  a flow  specification. 

Packets  matching  an  existing  classification  are  then  scheduled  to  continue  along  the 
flow  specification  with  the  associated  queue  prioritization  (Li  et  al.,  2000). 

INT-SERV  suffers  several  real-world  deficiencies.  It  requires  all  devices  in  the 
path  to  support  RSVP.  This  is  possible  to  achieve  in  a Local  Area  Network  (LAN),  but 
difficult  to  achieve  in  a public  Wide  Area  Network  (WAN).  RSVP  is  a half-duplex 
protocol.  This  requires  voice  conversations,  which  are  typically  bi-directional,  to  create  a 
flow  specification  for  each  direction  of  traffic.  INT-SERV  also  requires  a significant 
processing  effort  at  each  network  element  in  the  flow  specification  to  perform  the 
classification  and  scheduling.  If  it  were  to  be  used  on  major  ISP  trunk  connections, 
carrying  out  millions  of  packets  per  second,  the  overhead  per  packet  of  implementing  the 
necessary  checks  and  resource  management  is  almost  unacceptable  (Carpenter  and 
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Nichols,  2002).  INT-SERV  deficiencies  lead  to  the  generally  held  view  that  INT-SERV  is 
not  scalable  into  large  Voice-over-IP  networks  (Li  et  al.  2000).  As  summary,  INT-SERV 
is  not  realistic  in  high  demand  traffic.  It  is  not  possible  to  maximize  network  utility  while 
maintaining  certain  QoS  level  by  using  INT-SERV  since  it  wastes  capacity. 

2.3.2  Differentiated  Service  Model  (DIFF-SERV) 

Differentiated  Services  is  an  efficient  method  of  managing  traffic  based  on  its 
class.  DIFF-SERV  prioritizes  certain  types  of  traffic  like  voice  traffic  over  other  types  of 
communications.  It  works  by  categorizing  IP  packets  into  classes.  In  the  IPv4,  DIFF- 
SERV  uses  Type-of-Service  (TOS)  field  (one  byte  wide)  contained  in  the  IP  header  of 
each  packet  to  specify  a particular  behavior  type  that  determines  the  packet-forwarding 
scheme  and  priority.  In  IPv6,  the  traffic  class  byte  is  used.  DIFF-SERV  is  not  part  of  any 
negotiation  or  signaling  between  devices  themselves.  Priority  classes  are  assigned  by 
local  network  administrators.  The  assigned  tags  are  passed  in  the  packet  and  are  not 
subject  to  change  during  the  process  of  auto-negotiation  or  other  forms  of  signaling. 
They  help  administrators  prioritize  different  types  of  traffic  without  any  resource 
reservations  (Li  et  al.  2000  and  Transition  Networks,  2003). 

Differentiated  services  can  offer  the  following  classes: 

• Expedited  Forwarding  (EF),  which  defines  minimum  delay  and  jitter.  It  is  the 
preferred  mode  for  VoIP. 

• Assured  Forwarding  (AF),  which  introduces  three  selectable  packet  drop  rates. 
During  congestion,  packets  with  high  drop  precedence  are  discarded.  Thus 
enabling  the  more  important  traffic  marked  with  lower  drop  precedence  to  get 
through. 

• Best  effort  picks  up  the  remains  of  the  bandwidth  not  allocated  to  EF  and  AF. 
DIFF-SERV  requires  edge  routers  to  classify  traffic,  mark  the  traffic  for  rapid 

routing  decisions  by  core  routers,  and  police  connections  to  ensure  they  do  not  violate 
the  user’s  service  level  agreement  (SLA).  Core  routers  are  only  responsible  for  checking 
packet  precedence,  enabling  their  resources  to  be  directed  at  high-speed  switching. 
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Service  differentiation,  traffic  admission,  traffic  shaping,  and  other  computationally 
intense  functions  occur  at  the  network  edge.  DIFF-SERV  marks  the  TOS  field  based  one 
of  two  schemes:  quantitative  DIFF-SERV  or  priority-based  DIFF-SERV  (Li  et  al.  2000). 

In  quantitative  DIFF-SERV,  a deterministic  or  statistical  measure  of  throughput,  delay, 
jitter,  and  loss  are  calculated  to  provide  a quality  of  service  to  the  data  flow.  In  priority- 
based  DIFF-SERV,  data  flows  are  ranked  on  a relative  priority  scale. 

2.3.3  Asynchronous  Transfer  Mode  (ATM) 

ATM  offers  native  support  for  QoS  provisioning.  ATM  is  a very  successful 
protocol  in  multimedia  applications  because  it  can  provide  guaranteed  rates  and 
connections  that  are  so  valuable  to  real-time  traffic  such  as  voice  or  video  transmission. 
By  establishing  guaranteed  rates,  delay  and  jitter  can  be  reduced.  Since  ATM  provides 
guaranteed  rates,  packet  loss,  mostly  due  to  buffer  overflow,  is  a more  important 
parameter  than  average  packet  delay  (Kang  and  Tan,  1997). 

ATM  prioritizes  the  traffic  by  assigning  it  to  one  of  four  service  classes.  Each 
class  can  receive  a priority  level.  Four  ATM  service  priorities/queues  are  given  as  follows 
(Transition  Networks,  2003  and  The  ATM  Forum,  2002). 

• Constant  Bit  Rate  (CBR)  is  an  absolute  guarantee  of  a service  level.  This 
category  is  intended  for  real  time  applications  such  as  videoconferencing, 
interactive  audio  (e.g.  telephony),  audio/video  distribution,  (e.g.  television, 
distance  learning,  pay-per-view)  and  audio/video  retrieval  (e.g.  video-on-demand, 
audio  library). 

• Variable  Bit  Rate  (VBR)  is  for  variable  burstable  transmission  rates  with  a very 
good  throughput,  but  no  guarantee  as  to  the  consistency  over  time.  It  can  be 
used  for  streaming,  response  time  critical  transaction  processing  applications 
(e.g.  airline  reservations,  banking  transactions,  process  monitoring). 

• Available  Bit  Rate  (ABR)  offers  a minimum  guarantee  (Resource  exploitation, 
feedback  control,  critical  data  transfer). 

• Unspecified  Bit  Rate  (UBR)  makes  no  guarantees,  whatever  bandwidth  is  left  can 
be  used  (Best  effort).  Typical  applications  for  UBR  are  text/data/image  transfer, 
messaging,  distribution,  retrieval,  and  remote  terminal  (e.g.  telecommuting). 


21 


Three  multi-service  models  have  been  defined  above.  The  main  idea  among 
them  is  to  provide  different  treatment  for  different  classes  of  traffic.  These  methods  try  to 
maintain  appropriate  QoS  levels  by  using  priority-class  traffic.  Although  we  implement 
QoS  explicitly  in  our  model,  we  will  assume  that  the  underlying  network  handles  delay 
and  jitter  by  using  priority  classes. 


CHAPTER  3 
LITERATURE  REVIEW 

Existing  literature  on  QoS  can  be  classified  into  two  groups  for  the  purpose  of  our 
study.  The  first  group  considers  the  design  and  implementation  of  network  infrastructure 
and  operating  policies  including  QoS  issues.  The  main  concern  addressed  in  these 
problems  is  guaranteeing  a promised  level  of  quality  through  resource  allocation,  buffer 
management  and  scheduling.  The  second  group  in  the  literature  includes  pricing  and 
QoS  issues.  There  is  very  little  work  on  the  associated  business  models  and  the  problem 
the  customer  has  to  solve  to  optimize  its  position  within  the  context  of  a given  business 
model. 

There  is  extensive  work  on  QoS  in  electrical  engineering  and  computer  science 
fields.  We  have  already  covered  some  of  the  existing  literature  on  QoS  (Carpenter  and 
Nichols,  2002,  Habib  et  al.  2003,  Li  et  al.,  2000,  Miras,  2002,  Ragsdale  al.  2000, 
Rasmusson,  2002,  etc.)  that  describe  the  necessary  infrastructure  (including  buffer 
management  tools  and  packet  scheduling  algorithms)  and  service  models  (INT-SERV, 
DIFF-SERV  etc.)  that  support  applications  with  QoS  parameters.  Verma  et  al.  (1998) 
developed  a model  for  resource  reservation  at  nodes  along  the  path  between  source 
and  destination  node.  Their  solution  procedure  they  describe  is  similar  to  the  INT-SERV 
model.  After  deciding  which  call  is  admitted,  the  resource  is  allocated  throughout  path  by 
considering  end-to-end  delay  constraints.  They  describe  two  types  of  division  policies 
(reserve  of  transmission  rate  for  each  application)  to  allocate  the  resources  into  nodes. 
First,  equal  division  policy  where  the  receiving  node  decides  the  rate  for  each 
intermediate  node  transmission  rate  without  checking  any  load  on  that  node.  If 
intermediate  node  is  highly  utilized,  this  policy  may  result  in  higher  call  blocking  (blocking 
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of  new  call  attempts)  in  the  intermediate  node.  Therefore,  they  suggest  another  division 
policy  based  on  the  utilization  of  each  node.  They  show  that  their  algorithm  works  in 
heterogeneous  network  environments  with  multicast  traffic.  Yuang  and  Haung  (1998), 
Dou  and  Ou  (2001),  and  Pletka  et  al.  (2000)  studied  bandwidth  allocation  and  buffer 
management  for  DiFF-SERVE  models.  Yuang  and  Haung  (1998)  proposed  a bandwidth 
assignment  model  for  voice  and  data  traffic  to  guarantee  data  delay  and  voice  call 
blocking  probability.  They  used  a queuing  model  with  multiple  designated  channels  for 
data  and  voice  traffic.  They  argue  that  data  delay  increases  with  increasing  data  traffic 
intensity  and  mean  data  delay  is  highly  sensitive  to  the  ratio  of  voice  service  time  to  data 
service  time.  Dou  and  Ou  (2001)  proposed  five  bandwidth  reallocation  algorithms  for 
dynamic  provisioning  and  did  a performance  study.  They  argue  that,  when  most  of  the 
links  are  not  congested,  provisioning  should  not  be  restricted  to  congested  links, 
however,  for  the  congested  case,  provisioning  should  start  from  more  congested  links. 
Pletka  et  al.  (2000)  introduced  a new  buffer  management  scheme  that  can  support 
DIFF-SERV.  They  show  that  traffic  differentiation  is  possible  by  using  a delicate  packet- 
dropping mechanism  and  the  first  come  first  served  (FCFS)  dispatching  rule.  Their 
scheme  can  guarantee  a minimum  bandwidth,  support  excess  bandwidth  allocation,  and 
can  prevent  dropping  large  amounts  of  packets  from  a single  connection  by  using  drop 
precedence.  They  argue  that  with  this  scheme,  not  only  overall  buffer  usage  can  be 
optimized  packet  dropping  can  be  smoothed.  Low  and  Varaiya  (1993)  studied  service 
provisioning  in  ATM  Networks.  In  their  approach,  users  rent  resources  and  use  them 
based  on  their  own  traffic  and  quality  requirements  while  network  only  coordinates  the 
users’  choices  to  optimize  overall  network  performance.  However,  in  the  conventional 
approach  network  decides  resources  to  be  allocated  to  the  users.  They  consider  end-to- 
end  delay,  burstiness  of  messages,  and  different  type  of  service  requests  in  the 
provisioning  procedure.  Carpenter  and  Nichols  (2002)  claim  that  infrastructure  is  ready 
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to  use  for  QoS  supported  applications  and  more  products  with  some  level  of  DIFF-SERV 
support  are  entering  the  market.  Vendors  have  been  quick  to  incorporate  basic  DIFF- 
SERV  features  into  their  products.  The  usage  of  products  with  DIFF-SERV  functionality 
will  increase  in  the  near  future. 

Second  line  of  research  is  on  the  supplier  side  issues  including  pricing  and  QoS. 
In  this  line  of  research,  suppliers  are  assumed  rational  decision  makers  who  want  to 
maximize  their  utility.  Parameswaran  et  al.  (2001)  studied  a resource  allocation 
mechanism,  which  extend  the  DIFF-SERV  model  with  economic  functionality.  They 
showed  that  it  is  possible  to  combine  economics  and  network  engineering  to  derive  a 
feasible  economic  model  for  multi  service  architecture  in  spite  of  the  complex  nature  of 
the  traffic.  They  consider  allocation  and  auction  mechanisms.  In  their  formulation,  the 
market  (broker  company)  operates  the  mechanisms  to  compute  the  bandwidth 
allocations.  They  address  end-to-end  service  guarantees  with  constant  latency  and 
minimum  jitter.  They  discuss  that  differentiation  can  be  performed  by  using  user 
valuations  of  different  class  of  traffic  if  the  users  represent  their  need  for  service 
correctly.  Gupta  et  al.  (1999)  argue  that  a pricing  mechanism  can  be  used  as  a tool  to 
differentiate  different  traffic  requirements  considering  the  desired  QoS  of  customers. 
They  discuss  the  characteristics  of  feasible  pricing  mechanisms  that  support  multi- 
service traffic.  They  argue  that  (1)  pricing  should  not  be  at  the  packet  level  to  make 
implementation  easy  and  less  costly,  (2)  prices  help  congestion  control  by  encouraging 
users  use  network  when  it  is  less  congested,  and  (3)  multiple  priority  is  necessary  to 
take  care  of  different  applications  with  different  QoS  levels.  Reichl  et  al.  (2003)  reviewed 
the  Internet  pricing  models  and  proposed  a new  model  called  cumulus-pricing  model.  It 
is  based  on  a flat  rate  scheme  (fixed  charged  once  per  period,  usually  a month). 

Cumulus  pricing  model  specifies  types  of  services  with  expected  user  requirements  and 
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the  rate  required  to  get  each  service.  There  are  several  pricing  schemes  such  as  priority 
pricing,  real  time  pricing,  usage  based  pricing,  and  tax  band  pricing  etc. 

Einhorn  (1994),  and  Gupta  et  al.  (1997)  suggest  priority  pricing  as  one  of  the 
Internet  service  pricing  schemes.  Before  these  studies  priority  pricing  was  discussed  by 
Marchand  (1974)  for  queuing  systems  using  the  head-of-the-line  priority  rule  where 
customers  queue  in  the  system  according  to  priority  groups.  Einhorn  (1994)  mentions 
four  priority  classifications:  (1)  guaranteed  service  that  provides  certain  packet 
deliverability  or  guaranteed  delay  bounds,  (2)  predicted  service  that  provides  delay  and 
packet  loss  of  possible  magnitude,  (3)  link  sharing  service  that  allocate  bandwidth  to  a 
group  of  users  with  some  sharing  rules  and  (4)  datagram  service  which  has  no  quality 
commitment  except  best  effort  with  first  in  first  out  dispatching  rule.  He  argues  that  for 
the  network  providers  it  is  too  difficult  to  specify  in  advance  what  precise  menu  of  service 
offering  (charging  scheme)  is  appropriate  due  to  complexity  of  user  needs.  In  designing 
the  service  menu  each  provider  should  balance  the  benefits  and  costs  of  increasing 
service  diversification.  He  also  suggests  that  it  is  necessary  to  control  the  traffic  during 
some  specified  periods  such  as  busy  hours.  Users  should  be  able  to  reserve  bandwidth 
for  specified  periods  and  instantaneous  arrangements  should  be  available  for  short  run 
applications  as  well.  Service  offering  and  pricing  scheme  could  vary  by  bandwidth,  hours 
of  use,  and  delay  and  loss  probability  tolerance.  A user  is  interested  in  the  lowest  total 
cost  for  a service,  including  the  cost  of  delay,  since  it  maximizes  the  user’s  value  for  the 
requested  service.  Gupta  et  al.  (1997)  describe  the  priority  pricing  approach  for  modeling 
user  service  request  as  a stochastic  process  and  network  nodes  as  priority  queues.  A 
user  request  can  be  fulfilled  by  several  alternatives  (different  service  classes)  where 
each  alternative  has  different  price  and  expected  performance  (waiting  time  or  delay). 
Braun  and  Claffy  (1994)  and  Cocchi  et  al.  (1993)  discuss  that  it  is  possible  to  set  prices 
in  a network  with  multiple  service  priorities  so  that  users  of  every  application  type  are 
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more  satisfied  than  flat  pricing  considering  combined  cost  and  performance.  Therefore, 
priority  pricing  enhances  total  community  utility. 

Parris  and  Ferrari  (1992)  offer  a scheme  for  real  time  pricing  in  computer 
networks  that  allows  users  to  reserve  resources.  Charges  for  any  channel  depends  on 
resources  reserved,  type  of  service,  time  of  day,  and  channel  lifetime.  Type  of  service  is 
determined  by  considering  bandwidth,  buffers,  and  CPU  time  and  delay  resources 
reserved.  This  policy  is  unique  in  that  delay  is  characterized  as  a resource.  It  can  be 
considered  as  a quantification  of  the  loss  of  potential  for  the  network  to  admit  other 
channels  because  of  the  admission  of  the  channel  being  created.  Delay  resource  is 
inversely  proportional  with  delay  such  that  shorter  delay  requires  more  delay  resource. 
The  main  difficulty  for  real  time  pricing  is  that  it  is  difficult  to  reserve  resources  in  the 
current  Internet  (Braun  and  Claffy,  1994).  Mackie-Mason  and  Varian  (1993)  offer  another 
model  for  real  time  Internet  pricing.  They  argue  that  it  is  difficult  to  predict  in  advance  the 
optimal  price  of  network  service  due  to  internal  congestion  that  is  short-term,  bursty,  and 
unpredictable  phenomenon  in  most  part  of  the  Internet.  Their  proposed  pricing  scheme 
has  two  components:  (1)  a flat  connection  charge  that  is  based  on  the  type  of  customer, 
and  the  bandwidth  size  and  (2)  a congestion  charge  applied  during  network  congestion. 
They  claim  that  their  model  has  a long-term  advantage  that  it  properly  sets  congestion 
prices  for  valuing  capacity  expansion.  Efficient  pricing  of  network  congestion  in  the  short- 
term provides  the  optimal  investment  in  future  capacity.  According  to  Mackie-Mason  and 
Varian  (1994)  an  efficient  pricing  mechanism  should  have  a packet  charge  close  to  zero 
when  the  network  is  not  congested,  a positive  packet  charge  when  the  network  is 
congested,  and  a fixed  connection  charge  that  differs  from  user  to  user  (usage  sensitive 
charge). 

Altmann  and  Chu  (2001)  and  MacKie-Mason  and  Varian  (1995b)  studied  usage 
based  pricing.  They  argue  that  users  and  suppliers  prefer  a pricing  scheme  where  flat 
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rate  is  charged  for  basic  services,  usage  based  charge  is  applied  for  services  with  higher 
quality.  With  this  pricing  scheme,  users  will  benefit  both  using  unlimited  access  to  basic 
service  and  choosing  the  service  with  higher  QoS  when  they  need.  By  usage  based 
pricing,  congested  resources  can  be  prioritized  so  that  those  who  are  willing  to  pay  more 
get  the  higher  priority.  Therefore,  it  can  prevent  congestion.  Mackie-Mason  and  Varian 
(1995a)  have  found  that  the  equilibrium  price  and  capacity  will  maximize  net  social 
benefit  in  the  competitive  market  where  resources  is  provided  with  connection  fee  and 
usage  prices.  They  mention  that  due  to  externality  where  one  person’s  use  of  network 
resources  lowers  the  value  of  usage  for  everyone  else,  pricing  should  reflect  both  direct 
and  external  cost  of  usage  so  that  customers  will  use  the  resources  efficiently. 

Courcoubetis  et  al.  (2000)  have  studied  different  charging  schemes  with  effective 
bandwidth  as  basis  for  usage  charge.  Courcoubetis  et  al.  (1997)  argue  that 
determination  of  effective  bandwidth  can  be  problematic  since  it  is  hard  to  model  the 
moment  generating  functions  of  the  cell  arrival  process.  They  show  that  charging  for 
usage  with  effective  bandwidth  that  is  realized  bandwidth  for  the  customer,  can  be  well 
approximated  by  charges  based  on  simple  measurements.  Courcoubetis  et  al.  (2000) 
argue  two  services  can  be  substitutes  for  each  other  if  the  ratio  of  the  charges  for  these 
two  services  is  in  proportion  to  the  ratio  of  their  effective  bandwidths.  Effective 
bandwidth  is  a function  of  raw  bandwidth  and  the  statistical  characteristics  of  the 
connection  such  as  delay  and  packet  loss  ratio.  They  incorporate  delay  and  buffer 
overflow  (packet  loss)  QoS  constraints  into  their  formulation.  They  argue  that  network 
operator  will  provide  a small  set  of  QoS  differentiated  services  and  real  time  traffic  is 
offered  no  more  than  one  or  two  QoS  levels.  They  formulate  the  total  charge  as  a 
function  of  duration  of  connection  and  total  number  of  packets  carried.  They  also  discuss 
the  tax-band  charging  scheme.  They  argue  that  it  is  a suitable  scheme  in  case  of 
variable  rate  traffic,  so  that  a source  that  is  more  bursty  will  pay  a greater  charge.  Kelly 
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(1994)  also  formulated  the  pricing  structure  based  on  per  unit  time  and  cell  carried  over 
a call.  He  shows  that  effective  bandwidth  depends  on  user’s  expected  mean 
transmission  rate.  Courcoubetis  et  al.  (1996)  have  described  a pricing  scheme  for 
Available  Bit  Rate  Services  in  ATM.  In  their  model  prices  depend  on  demand  and  can  be 
iteratively  adjusted.  However,  in  best  effort  services,  prices  depend  on  average  cell  rate. 
According  to  their  scheme,  a connection  is  charged  by  considering  the  sum  of  the  price 
per  unit  of  bandwidth  on  all  links  along  its  path.  And  also,  users  are  free  to  adjust  their 
bandwidth  request  during  their  call  by  taking  into  account  unit  bandwidth  price 
announced  by  the  network  provider.  They  mention  that  in  the  case  of  optimal  social 
welfare,  the  prices  of  the  congested  links  maximize  social  welfare. 

Parris  et  al.  (1992b)  have  studied  different  pricing  schemes  including  peak  load 
pricing.  They  assume  that  QoS  is  negotiated  during  the  connection  setup.  The  network 
is  capable  of  making  QoS  guarantees  and  the  network  provider  charges  a price  for  each 
connection.  They  assume  QoS  requests  can  be  grouped  into  number  of  classes  such  as 
interactive  video  communication,  bulk  data  transfer  and  so  on.  They  also  mention  that 
the  price  can  be  made  sensitive  to  the  amount  of  reserved  resources,  the  number  of 
cells  or  packets  transmitted,  the  duration  of  connection,  the  time  of  day  and  the  quality  of 
the  service.  They  conclude  that  given  a pricing  scheme,  increasing  prices  first  increases 
then  decreases  the  net  revenue,  and  as  price  increases,  network  utilization  decreases. 
They  argue  that  peak  load  pricing  increases  revenue  by  spreading  demand  over  peak 
and  non-peak  periods.  Therefore,  it  is  useful  for  controlling  user  demand  to  prevent 
congestion. 

There  is  almost  no  work  about  the  customer’s  purchasing  strategy  when  faced 
with  a set  of  suppliers  potentially  offering  different  pricing  schemes.  There  are  some 
results  on  social  welfare  describing  when  a pricing  scheme  can  also  maximize  the 
customer’s  utility.  Thomas  et  al.  (2002)  show  the  existence  of  a social  welfare  (sum  of  all 
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user  benefits)  function.  They  determine  the  amount  of  required  resources  for  each  type 
of  service  in  order  to  maximize  the  sum  of  the  users’  utilities,  which  is  a reasonable 
description  for  a wide  variety  of  network  allocation  problems.  They  incorporate  QoS 
requirements  into  an  optimization  problem.  They  also  construct  an  algorithm  for  finding  a 
welfare  maximizing  solution.  Mackie-Mason  and  Varian  (1995a)  also  study  social  benefit 
maximization  as  described  earlier.  They  show  that  net  social  welfare  is  maximized  if  the 
resources  are  provided  with  connection  fees  and  usage  prices  in  the  competitive  market. 
The  usage  price  that  maximizes  profit  for  a monopoly  is  different  than  a social  optimum 
price. 

Though  not  directly  related  to  networking  there  is  some  work  in  the  supply-chain 
literature  analyzing  quality  of  service  issues  that  have  some  interesting  implications  for 
the  network  capacity  acquisition  problem  we  will  be  formulating.  For  example,  Chen  et 
al.  (2001)  have  studied  the  inspection  issue  for  production  inventory  systems.  They 
analyze  the  trade-off  between  buying  perfect  units  at  a higher  cost  and  buying  lower 
quality  units  while  spending  more  on  inspection  and  repair.  They  show  that  it  is  better  to 
inspect  a certain  number  of  units  and  repair  any  identified  defectives  rather  than 
inspecting  all  units.  As  for  what  quality  items  to  inspect  they  conclude  that  it  is  better  to 
inspect  the  units  with  a middle  quality  range.  They  argue  that  the  quality  of  units 
received  from  the  supply  source  is  improved  by  inspection  and  rework.  They  also 
suggest  a penalty  cost  for  unfulfilled  demand  for  each  unit  of  shortage.  In  our  model, 
there  is  an  analogy  between  loss  in  capacity  due  to  retransmission  and  rework  in 
manufacturing.  Especially,  inspection  is  an  interesting  problem  since  inspection  of 
quality  on  a network  requires  additional  resources  and  causes  further  delay.  The  same 
problem  of  how  frequently  to  inspect  and  what  item  classes  to  inspect  has  not  been 
discussed  in  the  telecommunications  literature  (since  most  traffic  to  date  has  been  best- 
effort). 
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Burke  et  al.  (2003)  study  the  supply  chain  sourcing  strategies.  They  focus  on 
supplier  selection  and  quantity  allocation  decisions.  They  characterize  conditions  under 
which  single-sourcing  and  multiple-sourcing  strategies  are  optimal.  Their  result  indicate 
that  multiple  sourcing  is  an  optimal  strategy  in  all  cases  except  when  supplier  capacities 
are  large  relative  to  the  product  demand  and  when  the  company  does  not  obtain 
diversification  profits.  In  this  case,  single  sourcing  is  the  dominant  strategy.  Since  there 
are  differences  between  the  supply  chain  and  the  telecom  network  environment,  very 
few  of  the  ideas  in  supply  chains  can  be  incorporated  into  our  model. 


CHAPTER  4 

VARIOUS  FORMULATIONS  BASED  ON  DIFFERENT  CHARGING  SCENARIOS 
We  formulate  a cost  minimization  problem  in  which  the  customer  strives  to 
acquire  network  capacity  subject  to  OoS,  capacity,  and  due-date  constraints.  We 
assume  that  an  environment  exists  in  which  the  customer  can  buy  network  capacity  with 
different  service  qualities  from  different  suppliers  who  can  offer  any  quality  of  service 
and  some  capacity  at  competitive  prices.  Therefore,  the  customer  has  access  to 
resources  with  different  QoS  levels,  capacities  and  pricing  schemes  and  has  tasks  to 
perform  with  different  QoS,  capacity  and  due-date  requirements. 

4.1  Definition  of  Service  Types  and  QoS  Requirements 
Generally  we  can  define  seven  different  service  types  for  different  applications. 
The  QoS  requirements  for  each  service  type  are  given  in  Table  4-1  that  is  an  extended 
version  of  the  Table  in  (Transition  Networks,  2003).  It  describes  how  sensitive  the 
service  type  is  to  a particular  type  of  QoS.  For  example,  batch  applications  are  relatively 
sensitive  to  delay. 

We  consider  different  types  of  traffic  with  different  bandwidth  requirements  and 
delay  tolerances  as  described  in  (Transition  Networks,  2003).  Three  parameters  (delay, 
jitter  and  packet  loss)  are  considered  as  relevant  QoS  definitions.  For  each  service  type 
different  levels  (low,  medium  or  high)  are  defined  for  each  parameter.  The  packet  loss  is 
defined  as  low  or  medium.  After  certain  threshold  value  (e.g.,  10%  for  VoIP  applications) 
the  packet  loss  is  assumed  as  high. 

The  duration  of  batch,  chat/online  customer  service  and  real  time  browsing/E- 
commerce  applications  depends  not  only  on  the  user  but  also  on  their  network’s 
throughput  levels.  However,  conversely,  the  duration  of  teleconferencing,  voice/video 
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and  video  streaming  depends  only  on  the  user.  The  data  flow  column  of  Table  4-1 
describes  the  direction  of  the  data  transfer  between  the  source  and  the  sink  excluding 
the  overhead  messages  such  as  confirmation.  For  example,  voice/video  and  video 
streaming  applications  are  one-way  real  time  since  at  the  logical  level  the  source  actively 
transmits  and  the  sink  always  receives.  Batch  transactions  are  also  one-way 
applications.  However,  VOIP  for  telephone  and  teleconferencing  are  duplex  type  of 
applications.  Both  sides  (receiver  and  sender)  can  be  active  at  the  same  time  during  the 
transmission.  In  the  browsing/E-commerce  applications  and  chat/online  customer 
service  applications,  the  interaction  mode  is  half  duplex  because  both  parties  usually 
wait  for  each  other’s  response. 

Table  4-1.  QoS  requirements  for  tasks 


Task  Type 

Throughput 

Latency 

Jitter 

Packet 

loss 

Impact 

Data 

Flow 

Task  Size/ 
Duration 

Batch/Data 

L 

L 

L 

T^ 

S“ 

F®7c 

Real  Time 
Browsing 
/E-Commerce 
Applications 

L 

L 

L/M 

T/PQ 

HD® 

cTJc 

Chat/Online 

Customer 

Service 

L 

H 

L 

L 

T/PQ 

HD 

F/C 

VOIP  (Phone) 

L 

H 

H 

L/M 

PQ'” 

D 

C/F 

VOIP 

(teleconference) 

M 

H 

H 

L/M 

PQ 

D 

C/F 

Voice/Video 

H ^ 

M/H 

H 

L/M 

PQ 

S 

C/F 

Video  Streaming 

M 

H 

H 

L/M 

T/PQ 

S 

C/F 

1:Low  6:  Fixed 

2;  Medium  7;  Compressible 

3:  High  8:  simplex 

4:  Time  9:  Half  Duplex 

5:  Realized  Quality  10;  Duplex 

Since  most  applications,  except  batch  transactions,  have  a human  involved  the 
realized  quality  of  the  application  depends  on  the  QoS  levels.  QoS  of  a network  affects 
the  duration  of  batch  applications.  Similarly,  due  to  buffering,  queuing  and 
retransmission,  the  transmission  time  could  change  in  browsing/E-commerce 
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applications,  chat/online  customer  service  and  video  streaming  that  directly  impacts 
realized  level  of  quality  by  the  human. 

On  addition  to  the  usual  classification  of  tasks  we  further  define  two  categories 
that  have  a profound  impact  on  the  decisions  to  be  made.  We  define  a service  to  be 
“time-fixed  (real-time,  size  compressible)"  if  its  size  can  be  changed  without  impacting 
the  successful  completion  of  the  task.  Most  audio/video  applications  are  size 
compressible  since  humans  cannot  detect  small  changes  in  the  amount  sent.  A service 
is  size-fixed  (time  compressible)  if  it  is  unaffected  by  the  delay  in  its  transmission  and  if 
all  bits  have  to  be  transmitted.  Most  data  applications  such  as  file  transfer,  batch 
database  transactions  are  of  this  category. 

Even  though  a more  detailed  breakdown  of  tasks  and  related  QoS  requirements 
is  possible  as  given  above,  the  two-task  type  classification  is  sufficient  for  our  purposes. 
In  our  formulation,  we  assume  that  a service  is  either  time-fixed  or  size-fixed  but  not 
both.  Note  that  in  certain  applications  the  amount  of  data  sent  might  affect  its  quality. 
This  however  is  not  part  of  the  QoS  provide  by  the  supplier,  yet  the  QoS  level  offered 
might  force  us  to  further  reduce  the  amount  of  data  sent  further  degrading  the  quality  of 
the  application. 

A data  application  has  fixed  size  but  one  can  tolerate  a long  delay  before 
completion  (high  packet  loss  ratio  and  delay  can  postpone  the  completion  of  the  task).  In 
general,  it  is  insensitive  to  throughput,  latency,  jitter  and  packet  loss.  Conversely,  a 
videoconference  is  highly  sensitive  to  QoS  levels.  Resources  with  low  QoS  impact  the 
realized  quality  of  the  task  (such  as  noise  in  speech,  frozen  frames  in  video,  etc). 
Moreover,  the  transmission  rate  might  need  to  be  lowered  because  of  bandwidth 
constraints  further  degrading  the  quality  at  the  receiving  node.  In  general  audio/video 
applications  have  a desired  target  transmission  rate  at  the  receiving  node  and  any 
deviation  from  it  costs  the  user  (opportunity  cost). 
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4.2  Problem  Notation 

A summary  of  notation  is  given  below. 

I,  J : The  ordered  index  set  of  resources  and  tasks  respectively. 

At,  As  : The  indexed  set  of  tasks  with  fixed  transmission  time  and  fixed  size  respectively 
where  J = and  A^V\A.,.=Q!> 

Qj  : The  worst  quality  level  tolerated  by  task  j.  qj  =q{5j,aj)  J^.,cr.  are  delay  and 
jitter  for  task  j respectively. 

Qi  : The  quality  level  guaranteed  by  resource  i.  Q.  - 0(<5',.,cr,.)  where  are  delay 
and  jitter  levels  of  resource  i respectively. 

a.  : Transmission  efficiency  that  is  calculated  as  one  minus  packet  loss  percentage 
of  resource  i. 

E;  : The  start  time  of  availability  of  resource  i. 

Li  : The  end  time  of  availability  of  resource  i. 

P,.  : The  bandwidth  of  resource  i. 

Tj  : Due  date  of  task  j. 

c-  : The  unit  cost  of  using  resource  i with  quality  Qj  based  on  total  number  of 
packets  successfully  carried. 

c'  : The  unit  cost  of  using  resource  i with  quality  Q;  based  on  total  transmission  time. 

cf  : The  unit  cost  of  using  resource  i with  quality  Qi  based  on  the  bandwidth  p 

purchased  from  suppliers. 

c°j  : The  opportunity  cost  of  missing  target  transmission  rate  for  task  j. 
t : Discretized  time,  t = 1,2,3,  ...T 

: Target  and  the  minimum  transmission  rate  of  task  j at  the  receiving  node. 
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: Estimated  scheduled  transmission  time  for  time-fixed  tasks  j e Aj. 

Xj  ; The  (fixed)  length  of  task  j eA^.  (Kbytes) 

The  decision  variables  used  are: 

Xj  : The  size  of  task  j e Aj.  (Kbytes) 

Vi  ; A binary  decision  variable  representing  the  purchase  of  resource,  i. 
rj  : The  transmission  rate  of  task  j.  (Kbytes/sec) 

yy  : A decision  variable  which  is  one  if  task  j is  assigned  to  resource  i,  zero 
otherwise. 

yyt  : A decision  variable  which  is  one  if  task  j is  active  (transmitting)  at  time  t on 
resource  i,  zero  otherwise, 
tj  : Start  time  of  task). 

4.3  Objective  Functions 

The  objective  function  is  composed  of  two  types  of  costs.  The  first  one, 
discussed  in  the  next  section,  is  the  telecommunication  cost  representing  different 
charging  schemes  of  suppliers.  The  second  is  the  potential  opportunity  costs  such  as  the 
cost  of  missing  deadlines  or  violating  due  dates,  the  cost  of  degradation  in  quality  (which 
is  more  relevant  for  real  time  streaming  applications),  and  the  QoS  switching  cost  when 
the  decision  maker  intentionally  downgrades  the  QoS  level  or  switch  to  another  supplier. 
4.3.1  Charging  Schemes 

We  model  a variety  of  charging  schemes  as  a function  of  the  amount  transmitted, 
the  transmission  rate,  bandwidth  and  QoS.  The  basic  functional  form  is  C = f^{Xj,rj,fi.) . 

We  assume  that  the  impact  of  QoS  on  the  price  is  implicitly  represented  in  function  f 
rather  than  an  explicit  parameter  of  the  function.  For  example,  for  a linear  functional  form 
the  unit  cost  used  in  one  function  might  be  higher  for  higher  QoS.  Under  this  general 
case  we  investigate  four  different  charging  schemes  commonly  studied  in  the  literature. 


36 


We  assume  that  the  supplier’s  pricing  is  rational  at  least  in  terms  of  QoS  and  bandwidth 
offered.  For  example,  every  supplier  charges  more  for  higher  QoS,  but  the  same  may 
not  hold  across  suppliers  (i.e.,  a supplier  might  sell  capacity  of  a higher  QoS  at  a lower 
price  when  compared  to  another  supplier). 

4.3. 1.1  Pay  per  byte 

C,  =fi(Xj)  = c-Xj 


The  total  charge  of  a task  is  a function  of  the  total  number  of  packets 
successfully  carried  (Altmann  and  Chu,  2001,  Gupta  et  al.  1997,  Parris  et  al.,  1992b). 
Since  the  bandwidth  is  not  part  of  the  cost  structure  we  assume  that  the  supplier  can 
offer  as  much  bandwidth  as  the  customer  demands  up  to  a published  maximum 
bandwidth  without  altering  the  cost  structure.  Such  a structure  should  encourage  the 
customer  to  use  as  much  bandwidth  as  possible.  Under  this  scheme  the  total  cost  of 
using  the  resources  is 

Totaicost  = X Z + Z Z ■ 

iel  jeA^  iel  jeAj- 


4.3.1. 2 Pay  per  byte  and  usage  time 

c'Ar. +c,:*Ar.r.Qr,.,y  €4 


C2=fiixj,rj)  = 


J J 

+ JgA^ 
r.a. 


The  total  charge  is  a function  of  the  duration  of  a connection  and  the  total 
number  of  packets  carried  (Altmann  and  Chu,  2001,  Kelly,  1994,  Courcoubetis  et  al., 
2000).  Such  a pricing  scheme  motivates  the  customer  to  reduce  the  duration  of  the 
connection  and  the  amount  of  data  sent.  For  time-fixed  applications  it  is  linearly 
increasing  in  r and  is  convex  decreasing  in  r for  size-fixed  applications.  The  total  cost  of 
using  resources  is 
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Total  Cost  = X Z Z 


/e/  y€/i7- 


le/  7€/4^. 


r 

yij<^i—  + C.Xj 


4.3. 1.3  Tax-band 


cjAty  + c^AtjV.a.,  jeA,  if  ^.  < y?,. 
c\  — + c-x  jeA^ 

r.ai 

c\At.  + At. r. a,,  j e if  /?;  > 

C,-'  ^ + c^x  JeA^  if  y?i  > A 

rjtti 


where  p,.  is  the  bandwidth  the  customer  desires  to  buy  and  use  (Courcoubetis  et 
al.2000). 

Note  that  the  total  transmission  rate  should  be  equal  to  the  bandwidth  purchased. 
For  time-fixed  tasks,  the  cost  function  is  piecewise  linear  as  given  in  Figure  4-1 . If 
customer  desires  /?,“  (or  ),  the  cost  function  will  be  curve  A (or  B)  in  Figure  4-1 . When 

increases,  the  slope  of  first  part  of  the  function  (between  0 and  ) decreases. 

However,  the  slope  in  the  second  part  increases  because  of  diseconomies  of  scale.  The 
suppliers  try  to  discourage  buying  unnecessary  amount  of  bandwidth  (Courcoubetis  et 
al.,  2000).  For  size-fixed  applications,  the  cost  function  is  hyperbolic  in  bandwidth  used. 


Total  Cost 


Bandwidth  used 


Figure  4-1.  Tax-band  charging  scheme 


4.3.1. 4 All-you-can-send 
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Q =/(/?,)- 


C: 


if 


In  this  case,  customer  will  buy  bandwidth  at  a fixed  price  (pay  per  resource)  for  a 
fixed  amount  of  time  (Altmann  and  Chu,  2001,  Gupta  et  al.,  1997,  Courcoubetis  et  al., 
2000).  This  pricing  structure  encourages  the  customer  to  use  as  much  bandwidth  as 
possible  in  order  to  maximize  his  utility.  The  cost  function  is  depicted  in  Figure  4-2. 

According  to  Figure  4-2,  if  the  customer  buys  a resource  with  bandwidth  , he 

has  to  pay  cf  regardless  of  how  much  of  the  bandwidth  he  uses.  Similarly,  one  has  to 

pay  cf  for  p\ . The  challenge  for  the  customer  is  to  acquire  just  enough  capacity. 


Total  Cost 


C] 


Bandwidth 


Figure  4-2.  All-you-can-send  type  charging  scheme 
4.3.2  Opportunity  Cost  Due  to  Quality  Loss 

A quick  inspection  of  the  objective  functions  and  the  constraint  sets  (below) 
shows  that  a formulation  with  only  acquisition  costs  encourages  sending  as  few  data  as 
possible  to  minimize  acquisition  costs.  In  reality,  firms  would  like  to  acquire  enough 
resources  while  maintaining  a level  of  quality  for  tasks  that  depend  on  data  sent  on  the 
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network.  To  incorporate  the  loss  due  to  realized  quality  we  incorporate  an  opportunity 
cost  in  our  model.  It  is  in  essence  similar  to  the  penalty  cost  of  Chen  et  al.  (2001 ). 

In  real  time  applications,  the  user  at  the  receiver  end  has  some  quality 
expectations.  Any  deviation  from  this  expectation  causes  dissatisfaction  with  the  quality 
of  the  task.  The  audio/visual  quality  of  real  time  applications  depends  on  their 
transmission  rate.  By  increasing  (or  decreasing)  transmission  rate  we  may  increase  (or 
decrease)  the  applications  quality.  For  each  real  time  task  we  assume  that  there  is  a 
well-defined  desired  transmission  rate,  observed  at  the  receiving  node  (r)^).  The 
opportunity  cost  for  each  task  is  calculated  as  follows. 

Opportunity  cost  = jc" 

4.4  Constraint  Sets 

We  consider  two  constraint  sets  that  will  be  combined  to  create  different  problem 
formulations.  The  first  set  contains  QoS,  capacity  and  due  date  constraints.  The  second 
set  contains  constraints  necessary  for  the  all-you-can-send  pricing  scheme  that  is 
modeled  in  Chapter  5. 

4.4.1  Constraint  Set  SI 


SI  contains  thirteen  equation  sets  as  follows. 


Z ^jy<j  + Z ^A^/jyij  - A (A-  - A ) V/  6 / 

je/i,  jeAj 

(4.1) 

IV 

m 

< 

m 

(4.2) 

^Ai-yi/j^  V/e/ 

(4.3) 

f'j^iyu -'^1'  V/e/ 

(4.4) 

Xj=Xj 

(4.5) 
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(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


Xj,rj,t.,fi.>0 

yy,yy,e{0,\) 

These  sets  of  equations  model  the  capacity  requirement,  QoS,  resource 
availability,  maximum  bandwidth  usage  and  scheduling  requirements  of  the  problem. 
Equation  set  (4.1)  models  the  capacity  requirements.  Equation  set  (4.2)  is  for  the  QoS 
requirement.  The  resource  has  to  have  a better  quality  level  than  what  is  demanded  by 
the  task.  For  example  if  a task  can  tolerate  a level  of  3 it  can  be  assigned  to  resources 
that  have  a quality  level  of  1,  2,  and  3.  Even  though  we  model  delay  and  jitter  as  levels 
we  explicitly  model  the  loss  probability.  It  is  assumed  that  transmission  rate  of  a task  is 
much  less  than  the  capacity  of  resources.  Equation  set  (4.3)  guarantees  that  the 
transmission  rate  of  a real  time  application  has  to  be  greater  than  the  minimum 
acceptable  receiving  rate  after  considering  the  loss  probability.  Equation  set  (4.4) 
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guarantees  that  the  transmission  rate  of  a real  time  application  has  to  be  less  than  the 
desired  target  transmission  rate  after  considering  the  loss  probability.  In  real  time 
applications  we  assume  that  lost  packets  are  not  retransmitted.  Therefore  the 
transmission  rate  after  packet  loss  has  to  be  greater  than  the  minimum  acceptable 
transmission.  Ideally  the  transmission  rate  should  be  equal  to  the  target  rate  at  the 
receiving  node.  Any  rate  above  this  limit  has  no  impact  on  the  quality  or  the  duration  of 
the  application  (since  most  devices  have  buffers  periods  of  higher  transmission  might 
alleviate  some  problems  but  we  assume  that  such  requirements  are  factored  in  while 
choosing  Equation  set  (4.5)  shows  that  size-fixed  tasks  have  predefined  task  size. 

Equation  set  (4.6)  shows  that  the  tasks  whose  service  type  in  the  time  fixed  set  have  a 
predefined  scheduled  transmission  time  that  is  constant  for  each  task.  Equations  set 
(4.7)  and  (4.8)  are  timing  constraints  and  indicate  that  each  task  has  to  be  allocated  to 
resources  when  they  are  available.  Equation  set  (4.9)  is  for  due  date  requirements  and 
indicates  that  each  task  has  to  be  allocated  to  a resource  and  completed  before  the 
deadline.  Equation  set  (4.10)  enforces  the  traffic  flow  at  any  given  time.  The  sum  of  the 
transmission  rate  of  all  active  tasks  on  each  resources  at  any  given  time  has  to  be  less 
than  or  equal  to  the  bandwidth  of  that  resource.  Equation  set  (4.1 1)  indicates  that  each 
task  will  be  assigned  to  only  one  resource.  Equation  set  (4.12)  and  (4.13)  guarantees 
that  the  jobs  are  actually  allocated  to  time  slices. 

The  term  g(Xj,r.,Q.)  in  equation  sets  (4.8)  and  (4.9)  is  the  expected  time  of 

transmission  for  task];  allocated  to  resource  i with  QoS  level  Qi.  It  can  be  formulated  as 
follows. 


j 


(4.14) 
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h{Si,a-)  is  the  probability  density  function  (pdf)  of  delay  Di  with  mean  5i  and  variance  Oi, 

(4.14)  describes  the  increase  in  the  completion  time  of  task  due  to  delay  and  jitter.  If  we 
only  consider  the  expected  value  of  Diwe  can  write  equation  (4.14)  as  follows. 


Unfortunately,  Equation  (4.15)  underestimates  the  actual  delay  of  a task  since 
most  often  for  a task  to  be  complete  all  its  packets  have  to  be  delivered.  Let  D.  be  the 

in  I 

delay  of  n‘*’  packet  of  task  j.  Then  the  delay  of  an  actual  task  over  resource  i is 

D.  - max„,  . Therefore,  the  expected  time  of  transmission  of  task  j will  be  given  as 

follows. 


Oj  is  an  index  representing  the  n'^  packet  of  task  j.  The  expected  time  of  transmission  of 
task]  increases  due  to  delay,  loss  probability  and  resubmission.  Since  VoIP  (both  phone 
and  teleconferencing)  and  VideoA/oice  service  types  are  real  time  applications  and  are 

not  size  fixed,  we  can  formulate  g[xj,r.,Q."j  as  follows. 


Note  that  this  does  not  impact  the  structure  of  (4.10)  since  ai  is  a coefficient 
rather  than  a decision  variable.  Delay  causes  some  interesting  complexity  in  modeling 
real  time  applications.  If  the  interaction  is  one-way,  delay  causes  a shift  in  the  broadcast. 


(4.15) 


X 

However,  if  the  interaction  is  two-way,  the  duration  of  the  task  is  -^  + KSi  where  K/2  is 


the  number  of  times  each  node  transmits  data  during  the  course  of  the  conversation 
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assuming  they  take  turns  in  transmitting  in  a Q&A  format.  We  will  not  pursue  explicit 
modeling  of  delay  and  jitter  in  the  scope  of  this  dissertation. 

4.4.2  Constraint  Set  S2 


S2  contains  two  equation  sets  as  follows. 


yij<v. 

(4.16) 

yij,  ^ yij  e I’  V/'  e y, r = E....L. 

(4.17) 

^i^ytj^yy,  ^ {0.1} 

Equation  set  (4.16)  shows  whether  resource  i is  being  used  or  not.  Equation  set 
(4.17)  shows  that  whether  at  least  one  task  has  been  assigned  to  resource  i or  not. 

4.5  Problem  Sets 

We  construct  four  different  problems  based  on  the  four  objective  functions 
discussed  in  section  4.3  and  constraint  sets  in  section  4.4.  To  simplify  the  analysis  we 
have  restricted  the  objective  functions  to  the  cost  of  acquisition  (telecommunication  cost) 
without  opportunity  costs.  In  Chapter  5,  we  will  focus  on  problem  C4S12  and  formulate 
the  problem  as  a trade-off  between  acquisition  and  opportunity  costs.  In  Chapter  7,  we 
will  focus  on  problem  C3S1  and  model  the  problem  as  a trade-off  between  acquisition 
and  opportunity  costs. 

4.5.1  Problem  C1S1 

Problem  Cl  SI  is  formulated  as  follows. 

min  X Z ^^jya  + Z Z 

le/  7'e/I,  iel  jeAj- 

subject  to  SI 

In  essence  problem  Cl  SI  is  similar  to  the  Generalized  Assignment  Problem 
(GAP)  and  is  NP-hard.  Ignoring  all  the  constraints  except  (4.1)  and  (4.1 1)  of  constraint 
set  SI  and  restricting  1=2  we  get  the  following  problem,  which  we  will  call  PO. 
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min  X Z ^^jya  + Z Z 

iel  JeA,  iel  jeAj 

subject  to 

J^x,y,.<fi,{L-E,)  '^iel 

jeJ 

Zyij=i> 

le/ 

e{0,l},  j3.,Xj>0 

Theorem  4-1:  Problem  PO  is  NP-hard. 

Proof: 

First,  we  can  observe  that  in  order  to  minimize  the  objective  function,  rj  cannot  be 
greater  than  rj' . Since  rj  is  the  only  decision  variable  on  the  second  part  of  the  objective 

function,  we  can  only  minimize  it,  if  we  choose  the  lowest  value  of  r,  that  is  rf . 

Second,  in  order  to  minimize  the  objective  function  we  need  to  maximize  'Y^Xjy.j 

J 

with  the  smallest  c/.  Therefore,  problem  Cl  SI  can  be  reduced  the  following  problem. 

max  ^ w-x- 

i 

subject  to 

i A(L,  - E,) 

JeJ 

6{0,1}  yjeJ 

This  is  a SUBSET-SUM  problem  and  it  is  known  to  be  NP-hard  (Martello  and  Toth, 

1990).  □ 
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4.5.2  Problem  C2S1 

Problem  C1S1  is  formulated  with  different  charging  scheme  (C2)  that  is  pay  per 
byte  and  usage  time.  The  objective  function  is  to  minimize  the  cost  of  using  resources. 
Problem  C2S1  is  formulated  as  follows. 


subject  to  SI 

4.5.3  Problem  C3S1 

In  this  case,  tax-band  charging  scheme  is  applied  to  problem  C1S1.  We 
postpone  the  explicit  formulation  to  Chapter  7 as  it  requires  additional  constraints  and 
variables  that  are  not  described  yet. 

4.5.4  Problem  C4S12 

Here  we  model  the  all  you  can  send  pricing  scheme.  Once  the  user  buys  the 
resource,  he  will  maximize  the  transmission  rate.  The  problem  C4S12  is  very  similar  to 
the  two-dimensional  bin-packing  problem  (2D-BPP)  and  formulated  as  follows. 


subject  to 


Vie/ 


(4.18) 


V/-€j 


(4.19) 


(4.20) 


rja^^y-jr/  \/j  e A^,  \fi&I 


(4.21) 


y/'ev4r,  V/e/ 


Xj=Xj  \/j^A, 


(4.22) 
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X; 


— = A?  . VJeAj. 


0 


yjeJ,  yiel 

yij^L.  VjeJ,yiGl, 
yy  ^ Tj  Vj  e J,  Vi  e / 

,/e^ 


tj+g 

1 

1 1 

tj^g 

yy<v.  VieiyjeJ 

yy,  ^ yy  e /,  vy  eJ,t  = E....L. 

Z Z yyt=^h’ 

iel  £,  Sf  <4 

Z H nyy<=^j’ 

/e/  Ei<t<Lj 


Xj,rj,tj>0 


e{0,l} 


(4.23) 


Bin  dimensions  are  raw  bandwidth,  p,.,  and  duration  of  resources  (Z,.  Note 

that  p,.  is  not  a decision  variable  in  this  formulation.  Tasks  are  represented  with  two- 

dimensional  shapes.  The  x-axis  and  y-axis  represents  usage  time  and  transmission  rate 
respectively  and  the  area  of  the  shapes  represents  the  size  of  the  tasks  (Xj).  Figure  4-3 
shows  an  example  of  bin  and  tasks. 

In  the  classical  bin-packing  problem  (Lodi  et  al.  2002),  the  size  of  the  items  and 
bins  are  fixed.  However  in  this  case,  the  x-dimension  of  items  representing  the  tasks  in 
set  At  is  fixed  but  the  areas  of  these  items  are  variable.  The  areas  of  the  items 
representing  the  tasks  in  the  set  As  are  fixed  but  the  shapes  of  the  items  are  variable. 
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Since  transmission  rate  for  time-fixed  tasks  is  variable  the  area  of  item  could  vary. 
Similarly,  the  transmission  rate  and  transmission  time  can  be  variable  unless  r.Ar^.  is 

fixed  for  the  size-fixed  tasks.  The  shape  of  items  in  the  set  As  could  be  a single 
rectangle  or  a combination  of  several  rectangles  such  as  task  number  5 in  Figure  4-3. 
Transmission  rate  of  size-fixed  tasks  can  change  during  transmission.  Constraint  sets 
(4.19)  and  (4.20)  satisfies  fixed  size  and  fixed  transmission  time  for  the  sets  As  and  Ay 
respectively.  The  y-axis  of  the  bins  is  controlled  by  the  constraint  sets  (4.21)  and  (4.22). 
The  x-axis  of  bins  is  controlled  by  constraint  set  (4.23).  At  any  time,  the  summation  of 
transmission  rates  of  active  tasks  has  to  be  less  than  or  equal  to  the  bandwidth  of 
resources.  Constraint  set  (4.18)  satisfies  the  summation  of  areas  of  items  has  to  be  less 
than  the  size  of  the  bins.  The  problem  reduces  to  first  scheduling  the  time  fixed  tasks 
into  resources.  Then,  allocated  the  size  fixed  tasks  into  available  spaces  in  the 
resources. 


Figure  4-3.  Bin  representation 

Theorem  4-2:  Problem  C4S12  is  NP-hard  in  the  strong  sense. 

Proof: 

Consider  a relaxation  of  the  constraint  set  S2  in  which  due  date  and  QoS 
constraints  are  relaxed  by  setting  Q=1  and  T-  =<x> . Assume  also  that  unit  cost  of  each 

resource  is  the  same  where  = 1 V/  e / and  all  tasks  are  time  fixed  (i.e.,  4=0). 
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At  optimality  r,.  = rj'  since  this  will  reduce  the  number  of  bins  required.  Consequently, 
we  can  rewrite  the  problem  as  follows. 

min^v,. 

(€/ 

subject  to 

Vi  el 

jeJ 

iel 


yy<v.  yiel,\/jej 

yy,  ^ yy  e A ^J,t  = ^,—A 

tj>E.yy  yjeJ,  V/€/ 


[a  +s(xj,  0 ’ Qi )]  yy  ^ A M/'  ^ e I 


Unyyi^^i  V/ e 7,r  = £,.... Z. 

jeJ 


^nyy^yy,^{^^) 

Finally,  the  two-dimensional  bin  packing  problem  is  known  to  be  NP-hard  in  the 
strong  sense  (Martello  and  Vigo,  1998;  Lodi  et  al.,  2002).  □ 

In  summary,  we  have  presented  the  notation,  the  charging  schemes  and  the 
constraint  sets  of  the  problem.  We  have  identified  4 possible  problem  types  considering 
combinations  of  constraint  sets  and  charging  schemes.  Interestingly  enough,  the 
problem  structure  changes  significantly  when  we  change  the  objective  function  or  when 
bandwidth  is  treated  as  a decision  variable  rather  than  data. 


CHAPTER  5 

CAPACITY  ACQUISITION  AND  TASK  ALLOCATION  SUBJECT  TO  QOS 
REQUIREMENTS  WITH  ALL  YOU  CAN  SEND  PRICING  SCHEME 


In  this  chapter,  we  formulate  a cost  minimization  problem  subject  to  QoS  and 
capacity  constraints.  For  tractability  we  relax  the  due  date  constraints  and  assume  that 
all  tasks  and  resources  are  available  at  time  zero  (Ej=0).  We  consider  the  all-you-can- 
send  pricing  scheme  in  which  the  customer  pays  a fixed  price  for  a fixed  bandwidth 
available  for  a fixed  duration.  We  also  assume  that  real-time  applications  have  desired 
transmission  rates  and  any  deviation  from  that  creates  an  opportunity  cost  for  the 
customer.  In  general  audio/video  applications  have  a desired  target  transmission  rate  at 
the  receiving  node  and  any  deviation  from  it  creates  some  dissatisfaction  for  the  user 
(opportunity  cost).  The  problem  that  the  customer  has  to  solve  is  a cost  minimization 
problem  that  reflects  a trade-off  between  cost  of  acquiring  resources  and  the  opportunity 
cost  of  degradation  in  realized  quality  that  is  important  in  the  real  time  applications. 


5.1  Problem  PI  Formulation 


(5.1) 


subject  to 


Z ^ jyij  + Z ^ ^ ^ 


(5.2) 


(5.3) 


(5.4) 


ZoTyv^/^i  V/e /,/  = !.. I, 


(5.5) 
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Eyij=i’ 

iel 

(5.6) 

V/-6  4 

(e/  tSLu 

(5.7) 

yy<Vi  \/ieiyjeJ 

(5.8) 

yyi^yy’  yieI,VjeJ,t^\..L. 

(5.9) 

(5.10) 

(5.11) 

rj,tj>0,  \/jeA^ 

(5.12) 

^nyy^yyr^{0A} 

(5.13) 

This  formulation  (P1)  considers  the  tradeoff  between  the  cost  of  acquiring 
capacity  and  the  cost  of  not  meeting  the  target  transmission  rates  in  real-time  tasks.  We 
assume  that  capacity  is  purchased  at  a fixed  price  for  a specific  bandwidth  and  duration 
(all-you-can-send  pricing). 

The  objective  function  contains  two  terms;  the  first  is  the  cost  of  acquiring 
capacity  and  the  second  is  an  opportunity  cost  due  to  not  meeting  the  target 
transmission  rate  in  real  time  applications  (i.e.,  time-fixed  tasks).  The  higher  the  capacity 
acquired  the  higher  is  the  acquisition  cost.  However,  purchasing  enough  capacity  allows 

for  setting  rjU.  = 7^.^  without  incurring  any  opportunity  cost.  Decreasing  the  transmission 
rate  of  the  tasks  increases  the  associated  opportunity  costs,  and  decreases  the 
acquisition  costs.  The  term  rjU.  represents  the  effective  transmission  rate  observed  at 
the  sink  node  after  taking  the  loss  probability  into  consideration.  Resources  with  low 
transmission  efficiency  (a.)  and  hence  lower  QoS  levels  usually  have  lower  prices  to  be 
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competitive  but  they  also  have  lower  effective  bandwidth.  Problem  PI  takes  this  fact  into 
consideration. 

Inequality  (5.4)  ensures  that  all  time-fixed  tasks  assigned  to  a resource  are 
completed  when  the  resource  is  available  (time  dimension).  Note  that  a similar  constraint 
is  not  needed  for  size-fixed  tasks  since  constraint  set  (5.2)  guarantees  that  a size-fixed 
task  is  only  assigned  if  there  is  enough  capacity  and  can  be  completed  on  time  since  we 
can  arbitrarily  vary  the  transmission  rate  (i.e.,  morph  its  shape). 

As  we  show  in  Chapter  4,  a special  case  of  this  problem  when  rj'  = r,.  = , 

cf  =c^ , Cj  =0  and  all  resources  are  identical  ( Q.  = Q,  L.  =L,  p.  = ^)  is  the  finite  two- 

dimensional  bin-packing  problem  known  to  be  NP-hard  in  the  strong  sense  (Martello  and 
Vigo,  1998;  Lodi  et  al.  2002).  Note  that,  if  all  tasks  are  size-fixed,  the  problem  turns  into 
one-dimensional  bin  packing  problem.  Since  some  of  the  decision  variables  {Vi,yy,yy, ) 

are  integer  (binary)  and  the  others  ( rj,  tj)  are  linear  and  equations  (5.1)  - (5.5),  (5.10), 

(5.1 1)  are  nonlinear,  PI  is  a mixed  integer  nonlinear  optimization  (MINLP)  problem.  As 
we  first  solve  this  problem  heuristically  we  leave  the  constraints  of  the  model  non-linear. 
Next,  we  present  a heuristic  and  computational  results. 

We  would  like  to  note  that  the  constraint  sets  described  do  not  guarantee  that  a 
task  is  completed  without  interruptions.  That  is,  once  a task  starts  its  transmission  rate 
might  drop  to  zero  for  dt  time  units  and  then  start  again  at  the  same  rate.  Even  though 
this  may  not  be  acceptable  in  real-life  it  is  a feasible  solution  in  this  formulation. 

It  might  be  possible  to  fix  the  solution  by  shifting  the  shapes  in  the  same  bin  to  eliminate 
the  interruption.  For  example,  consider  a solution  where  a task  is  split  into  two.  If  there  is 
a sub-rectangle  adjacent  to  one  of  the  pieces  while  including  the  other  then  it  is  possible 
to  fix  the  solution  (Figure  5-1)  by  swapping  the  sub  rectangle  and  the  piece  to  which  it  is 
adjacent.  However,  there  are  cases  where  a simple  swap  is  not  possible  (Figure  5-1). 
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Figure  5-1.  Swapping  of  tasks  to  eliminate  interruption  (a)  Solution  is  possible  by 
swapping  (b)  Solution  is  not  possible  by  swapping. 

Since  obtaining  an  exact  solution  is  relatively  hard,  simple  heuristics  can  be 
constructed.  We  discuss  one  such  heuristic  in  the  next  section. 

5.2  Heuristic  A 

Heuristic  A obtains  an  initial  solution  setting  the  transmission  rate  of  the  time- 
fixed  tasks  to  their  upper  bounds.  Under  these  conditions  the  opportunity  cost  is  zero 
but  the  capacity  acquisition  cost  is  maximum.  This  requires  using  the  largest  number  of 
resources.  The  main  idea  in  the  heuristic  is  to  reduce  the  cost  of  acquiring  capacity  by 
decreasing  the  transmission  rate  of  the  time  fixed  task  so  that  slack  capacity  could  be 
created  to  fit  some  other  tasks  (which  in  turn  increases  the  opportunity  cost).  When  it 
can  create  enough  slack  capacity  it  eliminates  one  of  the  used  resources. 

The  heuristic  is  described  below. 


Input:  vectors  T,  A/,  a,  q,  Q,  /3,  L,  c\  c“and  sets  A^,  I . Let  I.  be  the  set  of  jobs 
assigned  to  bin  i.  Let  B be  the  index  set  of  all  unused  bins,  C,.  = represent  the  total 
capacity  of  the  bin  and  L/,.  = C,.  - unused  capacity  of  bin  i. 

Output;  Total  cost,  bin  assignments  and  transmission  rates  (rj). 
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The  steps  of  heuristic  are  given  below. 

Step  1:  Calculate  effective  bandwidth, <- p^{\-a■),  available  bin 


C' 

sizeC,.  <-  C,  (l-c!r,  ),  and  unit  cost  cf  = — (assuming  the  bin  is  full)  for  each  bin  and 

Q 

sort  them  in  ascending  order  of  their  unit  cost  and  reorder  the  indexes  so  that  i=1  refers 
to  the  cheapest  bin. 

In  the  first  step,  the  unit  cost  of  each  bin  and  the  effective  bandwidth  of  each  bin 
are  calculated  by  considering  loss  probability. 

Step  2:  Find  the  available  bin  set  for  each  task  considering  QoS  requirements  of  tasks. 
Sort  them  in  decreasing  order  first  by  their  size  then  by  their  QoS  requirements. 

Step  3:  Let  ,j  g Aj.  . Assign  tasks  to  bins  starting  with  the  cheapest  bin  such  that 
A-  <-A-^{7}only  if  ^^  >0,. 


In  this  step,  time  fixed  tasks  has  been  scheduled  by  considering  their  upper 
bounds.  Heuristic  A schedules  the  task  into  the  first  available  bin  with  the  lowest  unit 
cost  if  the  QoS  level  of  the  bin  is  higher  or  equal  to  the  QoS  requirement  of  the  task. 


Step  4:  Compute  cf  = 


C: 


-,iel ! B . Sort  l/B  in  decreasing  order  of  cf  , Calculate 


the  current  cost  based  on  (5.1) 

In  this  step,  the  ratio  of  unit  cost  and  used  capacity  of  the  bin  are  updated. 
Basically,  if  a bin  is  almost  full,  cf  is  lower. 

Step  5:  Reduce  the  number  of  used  bins  if  it  is  feasible  and  economical.  If  the  solution 
changed,  go  to  step  4 otherwise  STQP.  Pseudo  code  of  step  5 is  given  below. 


i=1 


While  (i^H B){ 
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ye/,  ke//(Bi^i)  J^'^t 


[hen 


/<-/■  + 1 (look  at  another,  can’t  close) 


f \ 

Else  {Crash  jobs  in  At  such  that  '^Xj  < ^ Z^/  + Z rjAtj 

Je/i  ke/HB^i)  J^'^t  > 


(!-«.))  then 

j^Ai 


B<^B\Ji , compute  Ui^,c^  ,k  el  / B , 

Swap  tasks, 

sort  l/B  in  decreasing  order  of  l{Ci^-Ui^),keI ! B , i=1 

Else 


While  closing  the  bins,  we  consider  the  bin  with  a high  cf  ratio.  In  other  words, 
we  consider  the  bin  with  a high  unit  cost  and  low  usage.  Step  5 selects  the  bin  that  has 
the  highest  cf  value  (call  it  bin  1)  and  crashes  the  time-fixed  jobs  (reduces  the 
transmission  rate)  that  has  the  lowest  opportunity  cost  so  as  to  create  enough  slack 
capacity  to  fit  the  jobs  of  bin1  in  other  bins,  cf  is  designed  to  select  the  most  expensive 
bins  that  are  used  the  least.  For  example,  if  two  bins  have  approximately  the  same  per 
unit  cost  based  on  usage  (i.e.,  cf /(c,-t/,) ) Step  4 places  the  least  used  one  earlier  in  the 

order  as  given  in  the  following  example. 

Example  5-1: 


Bin  cost 


2000 


1000 
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Used  capacity  200 


^f/(C.-C/.)  10 

4l(C,-U,)'  0.06  0.1 

If  an  immediate  fit  is  not  possible  the  heuristic  removes  tasks  from  this  bin  until 
the  current  task  in  hand  can  fit  there.  This  can  be  carried  out  by  resolving  the  basic  bin- 
packing problem  exactly.  We  have  solved  it  using  a heuristic  that  swaps  the  jobs  among 
the  bins  so  as  to  fit  all  of  them.  In  order  to  remove  task  we  call  Swap  function.  The 
pseudo  code  of  swap  is  given  below. 

Let  be  the  index  set  of  all  unassigned  tasks. 

Sort  l/B  in  decreasing  order  of  Ui 
While 

For  each  y e 7 , do { 

For  each  / € 7/5  , do  { 
lf(7/,>x.) 

Assign  task  j to  bin  i,  I.  <- 1,  u {j},  <-  !{j} 

Sort  l/B  in  decreasing  order  of  U| 

Break 


Else 


while  ( Uf  < X-  and  {} ) { 

I;  <-  / {k},  <r-J^  u {k),x^  < x,,k,l  e /,.  and 
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(i.e.,  unassign  task  k starting  with  the  smallest 
until  there  is  enough  space).} 

/,.  <-  /,.  u {7},  <-  !{j}  (Assign  task  j) 

Sort  available  bin  size 
Break 

else  /<-/■  + 1 (go  to  the  next  bin) 

} 

} 

} 

} 


Swap  starts  with  the  largest  unassigned  task  and  the  bin  with  the  maximum 
availability.  If  an  immediate  fit  is  not  possible  the  heuristic  removes  tasks  from  this  bin 
until  the  current  task  in  hand  can  fit  there. 

Heuristic  A has  three  main  steps.  The  first  is  the  initial  assignment  of  tasks  to 
bins,  which  requires  solving  a 2D-BPP  (Martello  and  Vigo,  1998)  with  the  exception  that, 
jobs  in  As  can  take  any  shape  (r^.,  Ar^. ) as  long  as  the  area,  Xj,  remains  constant.  Note 


that,  any  solution  to  2D-BPP  is  a feasible  solution  to  PI.  Rather  than  optimally  solving 
this  initial  problem  we  use  the  rule  in  Step  3,  which  simply  picks  the  cheapest  bin  and 
assigns  as  many  jobs  as  possible  while  satisfying  the  QoS  requirements.  The  second  is 
the  rule  used  in  Step  5 to  select  a bin  to  close.  The  bin  that  has  the  highest  ratio  of 

Ck  ! [C,^  - V ,k  e I ! B\s  {he  most  expensive  and  least  used  bin.  This  makes  it  a 


perfect  candidate  for  closing  since  cf  gives  a good  measure  of  the  ratio  of  savings  to 
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opportunity  costs  due  to  closing  a bin  versus  crashing  jobs  (reducing  rj).  After  closing  the 

bin  with  highest  ratio  of  cf  /(Q  A:  e 7/5 , we  continue  try  closing  the  bin  with  the 

next  highest  ratio.  When  the  closing  of  the  bin  is  not  possible  for  a given  ratio,  the 
procedure  stop  because  next  available  bins  with  higher  ratios  also  be  impossible  to 
close.  The  third  important  step  is  the  crashing  rule  used  in  step  5 in  conjunction  with  the 

Swap  heuristic.  It  crashes  the  cheapest  tasks  (lowest  c° ) just  enough  so  as  to  fit  the 

unassigned  tasks  to  other  currently  used  bins. 

To  get  an  idea  of  the  quality  of  the  solution  we  also  compute  a lower  bound  (LB) 
that  is  given  in  the  next  section. 

5.3  A Lower  Bound  (Problem  P2) 

It  is  possible  to  obtain  a lower  bound  for  Problem  1 by  solving  the  following 
problem  (P2). 


(P2)  min  z = ^ cf 


jeAr 


iel 


iel  JeAj- 


(5.14) 


subject  to 


Z ^j-yu  + Z ^i^jya  ^ ^ ^ 


JqAj 


(q^-Q^)yy>0  yjeJ,  \fiel 

Zyij=^> 


/€/ 


0<yy<l,xj  <Xj 

where  Xj  = At j , x.  = r^-Atj , yij  represents  the  fraction  of  job  j assigned  to  bin  I 


(Figure  5-2),  Xj  is  the  amount  of  data  transmitted  and  c“  = — . The  objective  function 

At. 
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(5.14)  can  be  rewritten  as 


z = 


/ ye/lj'  i 


jeAr 


A simple  algorithm  (A1)  that  solves  P2  (ignoring  the  QoS  constraints  for  simplicity)  can 
be  specified  as  follows.Theorem  5-1  verifies  that  optimal  solution  of  P2  can  be  a lower 
bound  for  PI. 


A task  can  be  splitted 

/ 


Algorithm  A1: 

1 . Set  Xj  - and  sort  the  bins  in  ascending  order  of  cf' . 

2.  Assign  all  the  jobs  starting  with  the  first  bin  while  splitting  a task  if  it  does  not  fit  in 
the  bin. 

3.  Starting  with  the  job  that  has  the  lowest  c°  set  Xj  = x-  while  c"  < . 

Note  that  any  solution  to  PI  is  a feasible  solution  to  P2  since  P2  is  a relaxation  of  PI 
(the  integrality  constraints  for  yij  and  Vi  are  relaxed).  Consequently  the  solution  of  P2 
yields  a lower  bound  for  PI . 

Theorem  5-1:  A1  solves  P2 
Proof: 


Let  K be  the  index  set  of  the  selected  bins  when  the  algorithm  stops  and  let  K""  be 
its  complement.  By  definition  cf  <cf  ,i  e K,l  e K" . Let  also  M = {y|x^.  be  the  set 
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of  tasks  that  are  not  crashed  where  c®  > Max  cl . In  other  words,  M is  the  set  of  all  jobs 

whose  opportunity  costs  are  greater  than  the  maximum  acquisition  cost  of  the  bin  in  set 
K.  Let  also  N be  the  set  of  tasks  that  are  crashed  where  c°  < Max  cl . 

We  prove  that  the  sets  K,  M and  N yield  the  optimal  solution.  If  there  is  a better  solution 
it  can  be  obtained  either  by  using  a bin  that  is  not  in  K or  transferring  tasks  between  N 
and  M.  Since  K,  by  definition,  contains  the  cheapest  possible  bins,  swapping  one  unit  of 
capacity  cost  cj^  -cf  >0,k  e K,  I e . Therefore,  swapping  cannot  improve  the 
solution.  Consider,  uncrashing  the  most  expensive  task  in  N by  one  unit.  It  costs 
c/  -c"  > 0 to  take  this  action.  Similarly  crashing  a task  in  M costs  c°-cf  > 0 . 

Therefore,  further  crashing  or  uncrashing  tasks  cannot  improve  the  solution.  Hence,  A1 
solves  P2  optimally. 

In  the  next  section  we  present  computational  results  for  the  proposed  heuristic. 

5.3  Computational  Results  for  Heuristic  A 
We  generated  720  problems  with  36  different  problem  setting.  As  given  in  the 
Table  5-1 , we  have  identified  four  parameters  (pricing  structure,  tightness,  number  of 
tasks  per  bin  per  supplier,  and  ratio  of  the  number  of  time-fixed  tasks  to  number  of  size- 
fixed  tasks)  for  the  problem  setting.  We  have  randomized  all  other  factors  such  as 
minimum  and  maximum  bandwidth  requirement  for  time-fixed  tasks,  sizes  of  tasks,  bin 
capacities  and  QoS  requirements. 


Table  5-1.  Factors  in  problem  setting 


FACTOR 

LEVEL 

Pricing  Type 

3 

Tightness 

3 

# Tasks  Per  Bin 

2 

# TF  Tasks  / # SF  Tasks 

2 

Total  # Problem  Type 

36 
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Three  pricing  (one  random  and  two  supplier  sensitive)  scenarios  are  used.  In  the 
first  one,  all  bin  usage  costs  are  randomly  generated.  For  the  supplier  sensitive  cases 
(Figure  5-3),  the  bin  usage  cost  is  defined  as  a function  of  bandwidth  and  QoS  level  of 
the  bin.  It  has  been  generated  such  that  in  case  one,  there  is  a dominant  supplier  (i.e., 
the  cheapest).  In  case  two,  the  price  curves  of  different  suppliers  (as  a function  of  quality 
and  bandwidth)  intersect  each  other  so  that  at  some  bandwidth  and  QoS  level  a supplier 
becomes  more  expensive.  Intuitively,  the  customer  use  the  lowest  cost  as  long  as  the 
QoS  constraints  are  satisfied.  In  cases  where  QoS  cannot  be  satisfied  the  customer 
switches  to  a competitor.  Since  P1  is  a constraint  optimization  problem,  the  customer 
might  have  to  switch  supplier  due  to  QoS  requirements.  Most  economic  models  ignore 
such  constraints  and  may  not  see  the  same  behavior.  When  there  is  no  dominant 
supplier  though  the  customer  may  switch  to  a competitor  even  when  the  QoS  constraints 
are  not  violated.  This  should  occur  when  the  cost  of  using  additional  resources  from  the 
current  supplier  is  no  longer  the  cheapest.  We  study  the  random  cost  structure  solely  to 
test  the  performance  of  the  heuristic. 


Figure  5-3.  Supplier  sensitive  pricing,  (a)  parallel  cost  curves  (b)  intersecting  cost  curves 
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Another  factor  for  the  problem  is  tightness.  It  is  defined  as  the  ratio  of  the  total 
task  size  to  the  total  bin  capacity.  Three  levels  (50%,  70%,  and  90%)  of  tightness  have 
been  used  in  the  implementation.  It  is  expected  that  as  tightness  increases  the  number 
of  bins  used  and  the  number  of  crashed  tasks  increase 

The  two  levels  of  average  number  of  tasks  per  bin  arelO  and  25.  This  is  obtained 
by  setting  the  ratio  of  average  task  size  to  average  bin  capacity  to  0.1  and  0.04 
respectively.  We  randomized  the  sizes  of  tasks,  bin  capacities  and  the  QoS 
requirements  while  maintaining  the  above  ratios.  The  last  factor,  the  ratio  of  the  number 
of  time-fixed  tasks  to  that  of  the  size-fixed  tasks,  varies  from  1 to  2.  Note  that,  the  last 
three  factors  are  used  for  testing  the  performance  of  the  heuristic. 

Table  5-2  documents  the  average  solution  qualities  for  each  problem  setting.  The 
column  marked  Zpi/Zp2  represent  the  ratio  of  the  result  of  Heuristic  A to  the  lower  bound 
(LB).  The  overall  average  gap  between  solutions  of  Heuristic  A and  LB  is  4.9%.  The 
minimum  and  maximum  average  gap  is  2.4%  and  11%  respectively. 

The  maximum  gap  occurs  in  supplier  sensitive  pricing  with  intersecting  curve  as 
shown  in  Figure  5-4.  Because  of  the  relaxation  of  the  binary  constraints  on  y and  v in  LB, 
the  total  cost  curve  is  almost  linear.  Therefore,  the  gap  between  solution  in  Heuristic  A 
and  LB  may  sometimes  be  higher. 

Tables  5-3  and  5-4  and  Figure  5-5  show  the  effect  of  pricing  type.  When  the 
prices  are  random,  each  supplier  is  used  equally  (i.e.,  bin  usage  share  is  approximately 
25%  for  each).  However,  in  the  supplier  sensitive  pricing,  the  picture  is  quite  different. 
When  the  suppliers  have  parallel  cost  curves  (i.e.,  there  is  a low/high  cost  supplier  at  all 
levels  of  quality  and  bandwidth).  Supplier  1,  which  is  the  cheapest,  is  used  more.  Its  bin 
usage  cost  share  is  63.17%.  Intuitively,  we  need  to  stick  with  Supplier  1 and  its  bin 
usage  share  has  to  be  100%.  Since  this  is  a constrained  optimization  problem,  the 
optimal  solution  might  require  switching  suppliers  because  of  QoS  requirements.  For 
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example,  a single  supplier  may  not  satisfy  the  customer’s  quality  needs  even  if  it  is  the 
cheapest  supplier.  Another  reason  for  switching  is  due  to  the  combinatoric  nature  of  the 
problem.  Sometimes  using  a more  expensive  resource  instead  of  two  cheaper  resources 
might  reduce  wasted  capacity.  Most  models  discussed  in  Chapter  3 ignore  this  level  of 
detail  and  base  their  conclusions  only  on  price  or  quality. 

Table  5-2.  Average  solution  quality 


Zp2 

Zpi 

Zp2 

Zpi 

Problem 

Type 

Total 

Cost 

Total 

Cost 

Zpi/Zp2 

Problem 

Type 

Total 

Cost 

Total 

Cost 

Zpi/Zp2 

p1  t1  n1  r1 

273.75 

283.75 

1.0365 

p2  t2  n2  r1 

729.20 

756.40 

1.0373 

p1  t1  n1  r2 

281.74 

293.68 

1.0424 

p2  t2  n2  r2 

748.55 

783.94 

1.0473 

p1  t1  n2  r1 

254.32 

269.82 

1.0610 

p2  t3  n1  r1 

911.28 

955.59 

1.0486 

p1  t1  n2  r2 

271.09 

285.12 

1.0517 

p2  t3  n1  r2 

912.66 

970.42 

1.0633 

p1  t2  n1  r1 

388.37 

401.17 

1.0330 

p2  t3  n2  r1 

935.80 

967.29 

1.0336 

pi  t2  n1  r2 

405.69 

422.44 

1.0413 

p2  t3  n2  r2 

965.61 

998.09 

1.0336 

pi  t2  n2  r1 

377.88 

392.17 

1.0378 

p3  t1  n1  r1 

579.31 

612.83 

1.0579 

pi  t2  n2  r2 

391.31 

410.26 

1 .0484 

p3  t1  n1  r2 

595.53 

619.71 

1 .0406 

pi  t3  n1  r1 

494.12 

505.86 

1.0238 

p3  t1  n2  r1 

572.62 

601.43 

1.0503 

pi  t3  n1  r2 

528.72 

552.54 

1.0451 

p3  t1  n2  r2 

602.49 

641.64 

1.0650 

pi  t3  n2  r1 

493.24 

510.49 

1.0350 

p3  t2  n1  r1 

791.33 

831.59 

1.0509 

pi  t3  n2  r2 

522.52 

547.28 

1 .0474 

p3  t2  n1  r2 

807.93 

874.38 

1.0823 

p2  t1  n1  r1 

519.18 

534.17 

1.0289 

p3  t2  n2  r1 

838.80 

877.94 

1.0467 

p2_t1_n1  r2 

528.47 

544.14 

1.0296 

p3  t2  n2  r2 

855.58 

894.72 

1.0457 

p2  t1  n2  r1 

518.34 

542.38 

1.0464 

p3  t3  n1  r1 

1104.28 

1226.45 

1.1106 

p2  t1  n2  r2 

535.36 

564.05 

1.0536 

p3  t3  n1  r2 

1100.25 

1210.98 

1.1006 

p2  t2  n1  r1 

716.44 

743.32 

1.0375 

p3  t3  n2  r1 

1062.58 

1099.61 

1.0348 

p2  t2  n1  r2 

725.45 

764.26 

1.0535 

p3  t3  n2  r2 

1140.33 

1199.16 

1.0516 

Average  Zpi/Zp2 

1.0487 

Figure  5-4.  Maximum  average  gap  between  LB  and  Heuristic  A 


63 


Table  5-3.  Results  based  on  type  of  pricing. 


Pricing  Type 

Avg  # 
Bins 

Totai 

Cost 

# Crashed 
Tasks 

% Change  in 
the  Total 
Task  Size 

% Bin 
Utilization 

Run  Time 

Random 

7.75 

406.22 

35.54 

12.06 

98.96 

348.09 

Supplier  Sensitive  -A 

6.96 

760.34 

68.18 

18.51 

99.27 

158.77 

Supplier  Sensitive  -B 

7.22 

890.87 

75.22 

19.73 

99.18 

123.28 

In  the  third  case  where  suppliers  have  intersecting  cost  curves,  supplier  4 has 
more  bin  usage  share  as  85.73%.  Even  though  Supplier  1 is  cheaper  in  the  case  of  low 
bandwidth,  supplier  4 is  cheaper  for  high  bandwidth  and  high  quality. 

Table  5-4.  Bin  usage  cost  distribution  between  suppliers  considering  pricing  type. 


Pricing  Type 

Avg  # 
Bins 

Total 

Cost 

Supplier  1 
# Bin  Cost 

Supplier  2 
# Bin  Cost 

Supplier  3 
# Bin  Cost 

Supplier  4 
# Bin  Cost 

Random 

7.75 

406.22 

1.92 

90.48 

1.92 

90.13 

2.01 

94.43 

1.9 

88.53 

Supplier 
Sensitive  -A 

6.96 

760.34 

4.61 

426.75 

1.91 

198.04 

0.41 

47.58 

0.03 

3.19 

Supplier 
Sensitive  -B 

7.23 

890.87 

0 

0 

0.03 

8.95 

0.57 

104.32 

6.63 

680.76 

Figure  5-5.  The  effect  of  pricing  type 

The  results  of  the  effect  of  tightness  are  given  in  Table  5-5  and  Figure  5-6.  When 
tightness  increases  from  50%  to  90%  the  number  of  bins  used  and  the  total  cost 
increase  by  69.09%  and  88.61%  respectively.  Since  the  number  of  crashed  tasks  and 
the  percent  change  in  the  total  task  size  increase  as  tightness  increases,  percent  bin 
utilization  is  almost  the  same  in  all  tightness  levels. 
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Table  5-5.  Results  based  on  tightness 


Tightness 

(%) 

# 

Bins 

Totai 

Cost 

# Crashed 
Tasks 

% Change  in 
the  Totai  Task 
Size 

% Bin 
Utilization 

Run 

Time 

(Sec) 

50 

5.37 

482.73 

41.27 

11.69 

99.14 

231.99 

70 

7.47 

688.11 

61.26 

16.65 

99.22 

226.40 

90 

9.08 

910.50 

79.13 

19.85 

99.06 

187.27 

(a)  (b) 

Figure  5-6.  The  effect  of  tightness  on,  (a)  total  cost  (b)  % change  on  total  task  size 


Table  5-6  and  Figure  5-7  represent  the  effect  of  number  of  tasks  per  bin.  It  has 
almost  no  effect  on  number  of  bins  and  total  cost.  However,  there  is  less  of  a need  to 
crash  tasks  when  the  number  of  tasks  per  bin  increases,  since  the  size  of  each  task 
decreases  as  the  number  of  tasks  per  bin  increases.  Another  effect  of  smaller  task  sizes 
is  observed  in  the  utilization  column. 

Table  5-6.  Results  based  on  number  of  tasks  per  bin 


# Tasks 
Per  Bin 

# 

Bin 

s 

Total 

Cost 

# 

Crashed 

Tasks 

% Change  in  the 
Total  Task  Size 

% Bin 
Utilization 

Run  Time 
(Sec) 

10 

7.29 

685.96 

42.43 

19.14 

98.77 

31.18 

25 

7.34 

685.65 

76.86 

14.28 

99.48 

388.91 

Table  5-7  and  Figure  5-8  document  the  effect  of  the  ratio  of  the  number  of  time- 
fixed  tasks  to  size-fixed  tasks.  When  the  ratio  increases  from  1 to  2,  the  total  cost 
increases  3.84%.  This  is  due  to  the  fact  that  while  the  opportunity  cost  increases  (due  to 
more  crashing)  the  percent  bin  utilization  and  the  cost  remains  almost  the  same. 
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Figure  5-7.  The  effect  of  number  of  tasks  per  bin 


Table  5-7.  Results  based  on  the  ratio  of  the  number  of  time-fixed  to  size-fixed  tasks. 


TF/SF 

# Bins 

Total 

Cost 

# 

Crashed 

Tasks 

% Change  in 
the  Total  Task 
Size 

% Bin 
Utilization 

Run  Time 
(Sec) 

1 

7.30 

672.90 

50.31 

14.30 

99.19 

155.31 

2 

7.33 

698.71 

68.98 

19.08 

99.07 

264.78 

Figure  5-8.  The  effect  of  ratio  of  TF  / SF 

In  summary,  Heuristic  A seems  to  work  well  for  problems  with  random  pricing  but 
we  cannot  determine  its  quality  for  the  case  where  different  suppliers  are  cheaper  at 
different  bandwidth  ranges.  We  also  see  that  the  LP  relaxation  (P2)  is  tight  when  the 
resource  acquisition  prices  are  random  but  it  seems  to  degrade  when  there  are  suppliers 
dominating  at  different  bandwidth  ranges.  We  develop  a tighter  lower  bound  and  suggest 
a solution  procedure  in  the  next  chapter  and  test  Heuristic  A. 


CHAPTER  6 

A LOWER  BOUND  USING  THE  GENERALIZED  BENDERS  DECOMPOSITION 

METHOD 

In  the  previous  chapter,  we  formulated  the  cost  minimization  problem  subject  to 
QoS  and  capacity  constraints  and  solved  it  heuristically.  The  gap  between  the  lower 
bound  and  Heuristic  A might  either  be  due  to  a loose  lower  bound  or  a sub-par 
performance  by  Heuristic  A.  Hence  we  would  like  to  improve  the  lower  bound  further  to 
get  a better  feel  for  the  quality  of  Heuristic  A.  In  this  chapter,  we  develop  a new  lower 
bound  using  the  generalized  benders  decomposition  method  (GBD)  given  by  Geoffrion 
(1972).  Firstly,  we  present  the  lower  bound  formulation.  Secondly,  existing  solution 
methods  for  mixed  integer  problems  are  discussed.  Finally,  the  Generalized  Benders 
Decomposition  method  is  reviewed,  and  the  specific  adaptation  for  the  model  and  the 
implementation  results  are  presented. 

6.1  Lower  Bound  Formulation  by  using  GBD 
To  find  a lower  bound  (LB)  for  PI , we  consider  only  capacity  and  QoS 
constraints  and  relax  the  time  and  bandwidth  availability  constraints.  The  corresponding 
problem  (P3)  is  given  below. 

(P3)  Max  z = -cfv.  - X T-c°yy  + ^ ^^iCjyy 

iel  \ JeAj  JeAj  ^ 


subject  to 

Z ^jyij  + Z ^ A Vf  e / 

(6.1) 

je/iv  j^Aj 

(6.2) 
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(6.3) 

Ty^v,.,  ^JeA^,yieI 

(6.4) 

(6.5) 

7j‘-yy<rjai 

(6.6) 

(6.7) 

> 0 , Vy  e J 

(6.8) 

The  specific  adaptation  of  GBD  for  P3  is  given  in  section  6.4  after 
discussing  existing  solution  methods  and  giving  an  overview  of  GBD  in  the 
following  sections. 

6.2  Existing  Soiution  Methods  for  Mixed  Integer  Probiems 

In  a Mixed  Integer  Optimization  Problem  (MIP),  some  but  not  all  of  the  decision 
variables  are  restricted  to  integers.  Solving  such  problems  may  require  far  more 
computing  time  than  the  same  problem  without  the  integer  constraints.  Formulation 
given  in  Section  5.1  is  a mixed  integer  nonlinear  optimization  problem  (MINLP)  since  it 
has  at  least  one  nonlinear  constraint.  There  are  a variety  of  methods  suggested  in 
literature  dealing  with  “bin-packing  like"  problems.  One  of  the  well-known  methods  for 
solving  optimization  problems  is  the  lagrangian  relaxation  (LR)  method  (Fisher,  1981). 
Fisher  (1981)  reviewed  the  LR  method  for  solving  integer  programming  problems.  The 
idea  behind  LR  is  to  relax  (dualize)  some  of  the  constraints  and  add  them  to  the 
objective  function  with  a multiplier.  The  lagrangian  form  of  the  problem  is  usually  both 
easy  to  solve  and  its  optimal  value  is  a lower  bound  for  the  original  problem.  This 
method  has  been  used  in  many  different  optimization  problems  such  as  generalized 
assignment,  scheduling,  location,  partitioning  and  routing  problems  etc. 
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Martello  and  Toth  (1990)  give  both  approximate  (heuristics)  and  exact  algorithms 
for  BPP.  They  describe  heuristics  such  as  next  fit,  best  fit,  and  first  fit  that  are  used  to 
solve  bin-packing  problems.  They  also  give  a brief  overview  of  exact  algorithms  that 
exist  in  the  literature.  Exact  algorithms  are  basically  based  on  a branching  strategy  such 
as  the  first  fit  decreasing,  or  the  best  fit  decreasing  rule.  Lodi  et  al.  (2002)  present  a 
survey  of  the  two-dimensional  bin-packing  problem  (2D-BPP)  and  solution  procedures. 
They  mention  several  approximation  algorithms  such  as  next  fit  decreasing  height,  first 
fit  decreasing  height,  bottom  left  (packs  the  current  item  in  the  lowest  possible  position, 
left  justified),  floor  ceiling  (both  packing  the  items  from  left  to  right  with  their  bottom 
edges  on  the  level  floor  and  pack  the  items  from  right  to  left  with  their  top  edges  touching 
the  level  ceiling)  algorithms.  They  also  discuss  some  meta-heuristics  such  as  genetic 
algorithms  and  tabu  search  to  solve  2D-BPP.  Martello  and  Vigo  (1998)  derive  lower 
bounds  for  2D-BPP.  They  propose  an  exact  algorithm  in  which  they  have  embedded  the 
lower  bound  calculations  into  the  branch  and  bound  algorithm.  This  algorithm  is  the  first 
exact  algorithm  proposed  for  2D-BPP. 

6.2.1  Decomposition  Methods 

Dantzig  and  Wolfe  (1960)  have  proposed  a decomposition  technique  to  solve 
linear  problems.  The  Dantzig-Wolfe  approach  decomposes  the  problem  into  a series  of 
sub-problems  that  are  tied  together  by  dual  values.  A master  problem  is  created  which 
only  concentrates  on  the  coupling  constraints,  and  the  sub-problems  are  solved 
individually.  As  a result  only  a series  of  smaller  problems  need  to  be  solved. 

Benders  (1962)  generalized  the  ideas  of  Dantzig  and  Wolfe  to  solve  the  following 

MIP. 

min  cx  + dy 

subject  to 

Ax  + Dy>b 

x>0,  y>0,  X integer 
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Benders  decomposition  (BD)  algorithm  partitions  the  constraint  sets  and 
variables  into  two  exclusive  sets  and  divides  the  problem  into  two  parts,  a master 
problem  and  a sub-problem.  For  a fixed  value  of  the  integer  variables  the  sub-problem 
turns  into  a linear  problem  (LP).  After  solving  the  dual  of  the  LP,  the  dual  variables  are 
substituted  into  the  master  problem  and  solved.  At  the  each  iteration,  the  relaxation  gets 
tighter  and  the  lower  bound  increases.  BD  has  been  extensively  used  in  the  literature  to 
solve  mixed  integer  problems. 

Dantzig-Wolfe  method  separates  constraints  into  two  parts  (the  sub  problem  and 
the  master  problem).  In  contrast,  BD  separates  variables  into  two  parts.  Dantzig-Wolfe 
method  can  be  viewed  as  a dual  of  BD  or  vice  versa.  Applying  Dantzig-Wolfe  method  to 
the  primal  problem  is  equivalent  to  applying  Benders  to  the  dual  (Guo  et  al.  2002). 

6.3  Overview  of  the  Generalized  Benders  Decomposition 

Geofrion  (1972)  generalizes  BD  to  problems  where  the  constraint  set  does  not 
need  to  be  of  the  form  Ax  + Dy>b . We  call  this  procedure  GBD  and  first  give  a brief 
description. 

GBD  solves  nonlinear  programming  problems  (NLP)  of  the  following  form  . 

(NLP)  Maximize /(x, 

x,y 

subject  to  G(x,>;)>0,  (0.9)^ 

xeX,yeY 

where  y is  a vector  of  complicating  variables  because  of  the  fact  that  NLP  is  a much 
easier  optimization  problem  to  solve  in  x when  y is  temporarily  fixed  (Geoffrion,  1972). 
G(x,y)  in  (6.9)  is  an  m-vector  of  constraint  functions.  In  the  following  we  briefly  describe 


GBD. 
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Step  (1):  For  a given  y sY , solve  the  following  optimization  problem  that  is 
called  a sub-problem  (SUB)  and  obtain  a multiplier  vector  A associated  with  the 
constraints  G{x,y)  > 0 . Put  /['  = LB  = . 

(SUB)  Zsu3=  Maximize f{x,y) 

X 

subject  to  G(x,y)  > 0, 
xe  X 

Step  (2):  Solve  the  following  current  relaxed  master  problem  (RMP).  Let 
(j),  j)o)be  an  optimal  solution,  is  the  upper  bound  (UB)  on  the  value  of  the  optimal 
solution  to  the  problem  NLP. 

(RMP)  Maximize  y^ 

y^y,ya 

snh]QC\io y^<  Maximum  {f{x,y)  + A 'G{x,y)}J  = \X-  (6.10) 

x&X 

Step  (3):  Set  v = j)  and  solve  problem  SUB  again.  Set  LB  as  the  maximum  of  LB 
and  the  current  optimal  sub-problem  solution  (Zsue)  value.  W LB>UB-s , (where 
£■  > 0 is  a predefined  convergence  tolerance)  terminates  with  y and  the  corresponding 
solution  to  problem  SUB  as  the  solution  to  problem  NLP.  Otherwise,  go  to  step  (2)  with 
the  new  multiplier  vector  ZMhat  is  obtained  in  the  current  step).  If  revised  problem  SUB 
is  not  feasible,  optimal  multipliers  ) for  feasibility  cuts  can  be  generated  as  extreme 
rays.  In  this  case,  the  following  constraint  is  added  to  RMP. 

Maximum  { Ji^G{x,y)  } > 0,  j = 1,2,... 

xeX 

ZsuB  provides  a lower  bound  on  the  optimal  solution  to  problem  NLP.  At  each 
iteration,  problem  SUB  is  solved  by  using  different  y^^Y  and  LB  is  set  to  the  maximum 

of  these  sub-problem  solution  values.  The  optimal  multiplier  vector,  X , is  associated  with 
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the  constraints  G(;c,>^)  > 0 in  the  problem  SUB.  It  is  used  in  step  2 to  solve  the 
relaxation  of  the  main  problem  NLP. 

At  each  iteration,  the  problem  RMP  needs  to  be  solved  to  optimality  in  order  to 
determine  y to  be  used  in  step  (3)  for  finding  LB.  Note  that  a new  optimal  multiplier 

vector  is  found  at  each  iteration  j and  a new  cut  is  added  to  the  constraint  set  (6.10). 
The  upper  bound  is  updated  with  the  current  optimal  solution  y^  at  each  iteration.  In  this 

procedure  UB  is  nonincreasing,  since  each  problem  RMP  is  more  constrained  than  its 
predecessor.  Although  UB  in  step  (2)  is  nonincreasing,  Zsub  need  not  be  nondecreasing 
in  step  (3). 

Erenguc  et  al.  (1993)  have  used  GBD  to  solve  a project  scheduling  problem  that 
is  formulated  as  a MINLP.  Their  algorithm  performs  quite  well  even  when  the  problem 
size  is  moderately  large. 

6.4  Adaptation  of  the  Generalized  Benders  Decomposition  to  P3 

The  model  given  in  Chapter  5 is  a MINLP  with  complicating  variables  y and  v. 
Since  v is  dependent  on  y,  we  will  focus  on  y on  our  discussions.  Setting  _y  = 3; , the 

resulting  problem  is  linear.  Hence,  step  (3)  in  GBD  is  a LP.  However,  finding  y in  step 

(2)  is  a pure  0-1  Integer  Problem  (IP).  Before  applying  GBD  to  the  model,  it  is  better  to 
show  why  BD  would  not  work  for  this  model  in  the  next  section. 

6.4.1  Why  Benders  Decomposition  Method  does  not  Work  in  PI 
Since  PI  is  not  of  the  form 

maxlc'r  + /(y)jAr  + F(y)<b,  3;  e 

where  c and  b are  vectors  of  appropriate  dimension,  f is  a scalar  function  and  F(y)  is  a 
vector  function  of  appropriate  dimension,  BD  cannot  be  used.  The  reason  is  explained 


below  in  detail. 
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In  BD,  the  sub-problem  is  an  LP  for  a fixed  y . The  multiplier  set  U for  the  dual  of 
LP  has  to  be  independent  of>^ . In  this  case,  the  set  U is  dependent  on>^ , therefore,  the 
BD  method  is  not  applicable.  For  example,  let>»  = v and  consider  only  the  capacity 
constraint.  The  problem  reduces  to  the  following. 

Min  (6.11) 

(■  i jeA^ 

subject  to 

Z ^ jyij  + Z j^jyy  ^ A A e I 

ieA,  jeAr 


The  dual  of  the  problem  for  a fixed  y is  given  as  follows 


Max  2;  cfv.  + ^ ^ r,^'c°y,j  +2^/1,.  ^ x .yy  - a.p,L. 

ieJ  iel  jeAj  iel  jeA, 

subject  to 


Z^OA-Ty -Za°T!/'  V/'e4 

iel  iel 


As  seen  in  the  constraint  set  at  the  above  formulation,  the  multiplier  vector  Z is 
dependent  on  >^in  other  words  y . Since  A is  not  independent  oiy,  BD  method  is  not 
applicable  to  the  model. 

In  the  next  sections,  we  check  whether  problems  P1  and  P3  can  be  solved  by 
GBD  (Geoffrion,  1972). 

6.4.2  Checking  Property  P for  Problem  PI 

GBD  is  applicable  for  problems  where  Property  P holds  (Geoffrion,  1972). 
Property  P states  that  for  every  A > 0 the  supremum  of  / (a:,  y)  + A ^G{x,  y)  over  X in 

step  (2)  can  be  taken  essentially  independently  of  y.  For  simplicity  we  show  that  property 
P holds  on  a minor  relaxation  of  PI  (considering  all  packing  constraints),  but  the  time 
complexity  of  RMP  is  larger  compared  to  RMP  of  P3.  Consider  the  constraint  set  (5.2), 
(5.4)  and  (5.5)  given  in  Chapter  5. 
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le/ 


max  z = -cf  V,.  - r^'c^y^  + ^ 

J^At  jeAj- 


subject  to 

“Z  ^jyij  “ Z ^Ah'^iya  + ^ ^ 

ye  ^7- 


(6.12) 


^0  '^j&A^,y/i^i  (e.is) 

“Z%+A  V/e/,f  = 1..4  (6.14) 

j^J 

where  7 = {(v,,>;..,>;..,)|(5.3),(5.6),(5.7),(5.8),(5.9),(5.13)  Ao/r/}  for  fixed  r and  J 
Z = {(r^.,r.)  I (5.10), (5.11), (5.12)  hold] 

For  fixed  y eY  problem  is  a LP  in  xeX . 

Let  G{x,y)  > 0 s {(x,>^)  | (5.2), (5.4), (5.5)  hold] , f{x,y)  = objective  function  and 
X = vector  of  dual  variables.  Let  d.,i  e 7 be  dual  variables  corresponds  to  equation 
(6.12),  Uy,i  gIJ  eAj.  be  dual  variable  corresponds  to  equation  (6.13)  and 
w.,,/  el,te  L.  be  dual  variable  corresponds  to  equation  (6.14).  Then,  the  supremum  in 
step  (2)  can  be  written  as  f {y,X)  = SUP | / {x, >^)  + AG(x, >^)| . 


L'(y,X)^Y.  -cfv.  - 5]  r/c^yy 

iel  L JeAj- 


+ SUP< 

r,.l, 


Z Z wAij  - ^f>y>i  - (^^  • + ) y« + “//A- ) 


/€/  j^Aj 

+z^ 

iel 


■Z'^jyij-Z^A^jyij+^-M 


+ZZ^- 

iel  teL 


-Z  Gyij>  +A 

j^A-^ 
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r(y,^)-Z 


iGl 


'-cfy  + Z - liy  ) Yij  + “//'A  ]' 

JeA, 

+ SUP- 

iel  jeAj^ 

> 

+ z + Z ^./Ai 

^ J^'>s  '«A  / 

rjjj 

-^Ayyi-ZZ  Z 

le/  t^L,  jeAj- 

The  first  term  in  the  above  equation  is  not  dependent  on  r or  t.  For  the  second  term  we 
need  to  check  whether  Property  P holds  or  not.  In  order  to  maximize  the  second  term, 
tj  has  to  be  always  set  to  zero  since  -Uy  < 0.  In  that  case  the  solution  may  not  be 

feasible.  When  yy,  is  zero  for  all  t for  a given  i and  j,  yy  is  zero,  then  rj  becomes  a free 

variable  and  it  could  be  set  to  rj  = Tj’’ . However,  when  yy^  is  1 for  some  t,  yy  has  to  be 

1 . We  also  know  that  In  that  case  we  need  to  set  Vj as  given  below. 

I /si, 


For  y/  if  Max 


if  Max 


IJI 
/€i, 


—u 


f 


yt 

J 


> 0 =>  - Vj 


< 0 =>  r*  = rj' . 


We  give  brief  description  of  the  procedure  to  find  whether  q has  to  be  set  to  its 
upper  or  lower  bound  without  proof:  Fix  j,  for  each  i determine  the  Ar.  smallest 


A/, 


w.^  values  out  of  L|  such  values.  Let  w,.  = where  Wijq  means  the  t'*’  smallest  w^.  For 


/=! 


each  i compute  c°  -Atjdi~w..  If  the  maximum  of  these  is  positive  then  set  r,  to  its  upper 


A/, 

bound  else  to  its  lower  bound.  Computing  w.  for  each  bin  takes  0{L.  logZ,.)  (it 

t=\ 

is  possible  to  improve  this  further  though)  and  there  are  |l|  bins.  The  total  complexity  of 
this  step  is  o{^l\Max[L.  logZ.jj  which  is  pseudo-polynomial.  Consequently  we  will 
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focus  on  P3,  which  is  a simpler  version  of  P1 , in  the  hopes  of  developing  an  efficient 
exact  procedure  based  GBD. 

6.4.3  Checking  Property  P for  Problem  P3 

Let  X = [rj\{6.%)  hold}  7 = {(v.,;;..)  | (6.2), (6.3), (6.4), (6.7)  hold} . For  fixed 
yeY  problem  is  a LP  in  xeX . Let  G(x,>^)  > 0 = {(x,>^)  | (6.1),(6.5),(6.6)  hold] , 
f{x,y)  = objective  function  and  A = vector  of  dual  variables.  Let  u.,ie  I , 
cOy,i  G I,j  e Aj  and  ly-jA  e I,j  e Aj  be  the  dual  variables  corresponding  to  equations 
(6.1),  (6.5)  and  (6.6)  respectively.  Then,  the  constraint  in  step  (2)  of  GBD  can  be  written 
as  y^  <max|/(r,>;)  + /lG(r,>^)| . (Note  that  in  this  problem  a maximum  exists  so  we 

replace  the  supremum  with  maximum). 


Notice  that  some  of  the  terms  in  the  equation  (6.15)  are  independent  of  r..  These  terms 
are  constant  with  respect  to  Kj  . We  can  rewrite  the  equation  (6.15)  as  follows. 


To  ^ naax  + w.  ^ jYij  - X jrjYij  + A A 


(6.15) 


Z -Z  Z n [''/t,  -o«/] 


/ 


(6.16) 


In  order  to  solve  RMP  in  step  (2)  of  GBD,  first  we  need  to  maximize  r in  equation 
(6.16)  and  find  the  optimal  r*  value.  Theorem  6-1  describes  this  solution  showing  that 
property  P holds.  Notice  that  constraints  (6.5)  and  (6.6)  suggest  that  the  minimum  and 
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maximum  values  of  rj  depend  on  or,,  (i.e.,  rj'  < <7^^  if  assigned  to  bin  i). 

Consequently,  this  changes  the  exact  value  of  r* depending  on  the  bin  assignment. 
Fortunately,  since  the  job  can  be  assigned  to  only  one  bin,  the  or,,  in  Equation  (6.16) 
cancels  this  effect  which  enables  us  to  work  with  rl'  ,7^  values  independent  of  or,. . 

Theorem  6-1: 

The  solution  given  by 


solves 


if  Mu  (c“  - At jM.  - cOy  + y/.j ) < 0,  V/  =>  = r/ , Vy 

(6.17) 

if  M(^[c°  - At  jM;.  -ay+(p^j)>0,3i  =>  rJ  = 7^ , Vy 

(6.18) 

™ Z 0 Z [yy  ( ^ V/  -^y)  + ^y]\ 

/€/  -'J 

(6.19) 

Proof: 


Since  r,.  > 0,  or,.  > 0,  >0,  iel,j€j  , and  one  and  only  one  y.j  =1  for  a fixed 

j,  we  need  to  consider  the  cost  coefficients,  (c“  - At j-wj.  + , one  at  a time  based  on 

i (rather  than  the  sum  of  an  arbitrary  number  of  "i”s).  If  one  of  them  is  positive,  setting  rj 
to  its  maximum  maximizes  (6.19).  This  structure  enables  us  to  examine  only  a linear 
number  of  coefficients.  Since  each  rj  can  be  determined  independently  we  can  focus  on 
a rule  for  an  arbitrary  task  j.  In  short,  we  have  to  maximize  the  following  inner  summation 
term  for  each  y e ^ , 

0 Z (<  - (6.20) 

iel 

We  show  by  contradiction  that  rules  (6.17)  and  (6.18)  always  improve  the  optimal 
solution.  Assume  r,.  is  an  optimal  solution  where  0 < r,'  < 7^^ . 
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Let  {c°  - AtjM.  - co.j  +y/.j^>Q . If  we  set  r.  = , since  r-  < r^'  we  improve  the 

value  of  (6.20)  hence  the  solution  improves.  This  is  a contradiction  since  we  assumed 
that  Vj  is  optimal.  Let  (c°  - AtjM,.  < 0.  Equation  (6.20)  is  decreases  as  r 


increases.  Setting  -rj'  maximizes  (6.20)  creating  a contradiction. 


6.4.4  GBD  procedure  for  P3 


After  finding  the  value  r for  each  j,  we  can  plug  it  in  to  equation  (6.16)  again  and 
obtain  the  following  equation. 


/€/ 


+ Z - Z y^iiyv^j  - Z 

j&Ai  j^Af  j^Af 

Z X jyij  + Z “.T/,  -"//)+ 0 


-cfv. 


(6.21) 


V L J ^ 

Note  that  there  are  four  terms  inside  the  summation  in  equation  (6.21 ).  The  first  term  is 
constant,  the  second  term  is  linear  with  variable  v. , the  third  and  fourth  terms  are  also 


linear  with  variable  y.j.  The  third  and  fourth  terms  are  the  summations  over  As  and  Aj 

respectively.  Equation  (6.21)  is  the  constraint  that  is  calculated  and  added  to  RMP  in 
each  iteration  at  step  (2)  during  GBD.  By  adding  new  constraint,  we  are  adding  another 
cut  to  the  feasible  set  in  order  to  find  the  optimal  solution. 

6.4.4.1  Preprocessing  for  Step  (1) 

In  GBD,  initial  assignments  for  y and  v values  are  given  in  step  (1 ) and  sub- 
problem are  solved  depending  on  these  integer  value  assignments.  Since  the  dual 
multipliers  are  used  in  step  (2),  this  initial  assignment  might  be  important  for  the 
convergence  rate  of  GBD.  If  we  have  good  assignment,  we  may  have  the  optimum 
solution  for  the  master  problem  in  less  iteration.  For  the  initial  assignment,  we  are  using 
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the  same  initial  solution  used  by  Heuristic  A of  Chapter  5 with  the  exception  that  we  do 
not  have  some  of  the  constraints. 

G.4.4.2  Feasibility  Check  in  Step  (2) 

In  step  (2),  the  following  RMP  that  is  0-1  IP  is  solved. 

(RMP)  max 
subject  to 

i€l 

y^j<v.,  V/'G^,  Vfe/ 


qj>Q.,yy  VieiyjeJ 
Z + Z < a,f3,L.  V/  g / 


(6.22) 

(6.23) 


yo^ 


JgAj 

f 

u-a.p,L.  + 

S - 

V 

j^Aj 

i^Aj 

jeAr 

+ Z w 

At, 47* 

-1  \ 


-cfv. 


JsAt 


V/G/,y/'G  = 


^,=ne{o,i} 


Notice  that  there  are  two  types  of  solutions  from  the  problem  SUB  in  step  (3). 
Either  the  problem  has  a solution  that  is  feasible  or  the  problem  is  infeasible.  In  the 
former  case,  dual  variables  are  used  for  finding  optimality  cuts.  In  the  latter  case, 
extreme  rays  are  used  to  find  the  parameters  for  feasibility  cuts.  In  order  not  to  avoid 
extreme  rays,  we  need  to  guarantee  that  a solution  of  step  (2)  yields  a feasible  solution 
in  step  (3).  This  can  be  done  by  imposing  (6.22)  and  (6.23)  that  is  a modified  version  of 
6.6. 

6.4.4.3  Solving  Problem  SUB  in  Step  3 

In  step  (3),  the  following  SUB  is  solved. 
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r \ 

(P3)  max  z = ^ -cfv.  - c°yy  + ^ rja,c°yy 

iel  JeAj  jsAj  ) 

subject  to 

Z ^iya  + Z 6 ^ 

>€-4,.  ;e/)r 

V/'e4>  V/6/ 
r.yy<rya^  \/j  ^ A^,  \fi  e I 
rj>0  , \fj  e J 


Notice  that  when  v,.  = v,.,  i e I and  y^  = y^,  iel,j  eJ , SUB  becomes  an  LP.  In 

the  next  section,  the  GBD  method  is  illustrated  in  an  example. 

G.4.4.4  Example 


We  consider  3 resources  and  4 tasks  with  the  parameters  given  in  Table  6-1  and 
Table  6-2. 

Table  6-1.  Bin  parameters. 


* 

Qi 

Li 

Pi 

CiP 

OCi 

1 

2 

25 

35 

30 

0.99 

2 

1 

35 

20 

35 

0.98 

3 

1 

40 

30 

40 

0.97 

Table  6-2.  Task  parameters. 


J 

qj 

Type 

Atj 

A 

1 

1 

TF 

20 

15 

— 

20 

2 

2 

2 

TF 

15 

15 

- 

20 

1 

3 

1 

TF 

25 

10 

- 

15 

3 

4 

2 

SF 

- 

- 

300 

— 

0 
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The  sizes  of  the  bins  are  calculated  by  C,.  = PiL.y.  and  found  as  Ci  =866.25,  C2  =686 
and  C3  =1164.  We  formulate  the  problem  as  a mathematical  program  of  the  form  given 
by  PI. 


3 3 


3 3 


max  z = - J]  cfv.  - ^ Z + Z Z 


1=1 


<•=1  ;=i 


1=1  y=l 


z = -105-30v,  -35v2  -4OV3  +2r,  +«3>^3,) 

+ ''2  («iT|2  + «2^22  + «3T32  ) + 3^3  («lTl3  + ^2^23  + «3T33  ) 

subject  to 

(20r,j2,,  +15r2>^,2 + 25r3>^,3)a,  +300>’,4  <866.25 
(20r,;23,  +15r2j233  +25r3>^23)^i  +300^24  ^686 
(20r,>^3,  +15r2>^32  + 25r3>^33)cir,  +300^,4  <1164 

^ 22T21 

^2  ^ QT,2 
<h  ^ Qiyii 
^2  - Q^yii 

^3  ^ QT,3 
^3  ^ 02^23 
^3  ^ 03^33 

Qa  ^ Q,yu 

^4  ^ 02^24 
^4  ^ Q^yiA 

TU+T21+T31  =1 

T12+T22+T32  =1 

T.3+T23+T33  =1 
TM+T24+T34  =1 

Tu  ^v, 

T,2 

T,3  ^ 

Ti4  ^ 
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>^2.  ^ ^2 
>^22  ^ ^^2 
>^23  ^ ^^2 
>^24  ^ ^2 

>^32  ^ ^3 
>^33  ^ ^3 

;^34  ^ Vj 


ax^xyu 

«2''i3^2. 

«3'1>^31  ^'•'x" 

c^x^yxi  ^^2 
(^iriyii  ^^2 

a^r^y^j 

Uxr.yx,  < 7,^ 

(^2hy22^'T2 

C(2ryy,,<7,^ 


^X^X 

a^rx  >^^21 

a,r, 

«l''2  ^^Vl2 
*^2^2  - ^2  ^22 
«3''2^^^32 
a,r3  >^V,3 
^2^3  — ^3  ^23 
*^3^3  — ^3  ^33 

r.  >0 

STEP  (1):  we  begin  with  selecting  y eY  then  solving  the  initial  problem.  After 
preprocessing  y^^  = jVj,  = y^^  = y^4  = 1,  v,  = Vj  = 1 and  all  other  yy  and  Vi’s  are  set  to  zero. 
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Therefore,  we  solve  the  following  LP  problem  in  step  (1). 

max  z = -175  + 1.94r, +0.99^2 +2.91^3 

subject  to 

14.85^2  <566.25 
19.4r, +24.25^3  <1164 

0.99r2  < 20 
0.97r,  < 20 
0.97r3  <15 
0.99^2  >15 
0.97r,  >15 
0.97r3  >10 

After  solving  the  problem  we  obtain  LB  = -70,  ry,2  = 1,  = 2,  co-^^  = 3 and  all  other 

dual  variables  as  zero. 

ITERATION  1 

STEP  (2):  This  step  has  two  parts. 

Part  a:  First,  we  need  to  find  r*  values  considering  Theorem  6-1.  Since 
[c°-  AtjW^.  - w..  + ) > 0, 3/  for  all  tasks,  we  set  r*  = for  all  time  fixed  tasks.  In  other 

word,  we  set  r,’  = 20,  = 20,  r*  -\5 

Part  b:  After  finding  r , we  find  the  cut  and  solve  the  following  IP  problem. 

Max 

subject  to 

To  ^(40;^n  +40>^2.  +20>^22  +20>^32  +45>^,3  +45>^23  -30v,  -35v2  -4OV3) 

Xyi,=i.  v,sj 

iel 

Qj^Qiyy  ^ieiyj^J 
yy<v.,  VjeAj.,  Vi  el 

Z ^jyij  + Z < a.p,L^  Vi  e I 

jeAj- 
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After  solving  the  problem,  we  obtain  UB  - 20,  v,  = = 1,  y^2  = ~ = >^21  = ^ 

all  other  y and  v values  as  zero. 

STEP  (3):  In  this  step  we  substitute  the  v and  y values  found  in  step  2 in  the  main 

problem  and  solve  it.  The  formulation  of  LP  problem  is  given  as  follows. 

Max  z = -170  + 1. 96r, +0.992-2 +2.972-3 
subject  to 

14.85r2 +24.752-3  <566.25 

19.62-,  < 686 

0.992-2  ^ 20 

0.982-,  < 20 

0.992-3  <15 

0.992-2  >15 

0.982-,  >15 

0.992-3  >10 

After  solving  the  problem  we  obtain  LB  = -74.05,  w,  = 0. 12,  co^y  = 2,  = 0.8  and  all 

other  dual  variables  as  zero. 

ITERATION  2 

STEP  (2):  First,  we  set  r,’  = 20,  rl  = 20,  2-3*  = 15  in  part  a,  then  we  solve  the  following 
IP  problem  in  part  b. 

Max  >'0 
subject  to 

yo  ^ (40>'„  + 40>^2i  + 20>'22  + 20^32  + 45>',3  + 45>^23  - 30v,  - 35v2  -4OV3) 

^0  +40>^3,  -87,,  +20^32  +20>^22  +45^33  -28>^,2  +45>^23  -30v,  -35v2  -4OV3  +54.95) 

Xyij=l,  V/'eJ 

iel 

Qj^Qtyij  ^iel,yjej 

Vye^,  V/e/ 

Z + Z ^ V/  6 / 

je/(,  JSA, 

After  solving  the  problem,  we  obtain 

UB  - -29.05,  V,  = V2  = 1,  yy2  = yy^  = >^21  = >^23  = ^ all  other  y and  v values  as  zero. 
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STEP  (3):  When  we  substitute  the  v and  y values  found  in  step  2 in  the  main  problem 
and  solve  the  LP  problem,  we  obtain  LB  - -13.9,  =0.1,  cOj^  =0.5,  £0^2  =land  all 

other  dual  variables  as  zero. 

Note  that  UB  is  nonincreasing,  however,  LB  is  not  nondecreasing  (it  need  not  be 
nondecreasing).  We  found  the  optimum  solution  after  6 iterations.  At  the  optimum,  all 
tasks  are  assigned  to  the  third  bin  and  the  solution  for  the  main  problem 
A-*  = (20.6186,20.202,15.4639)  with  an  objective  value  of -70.  Note  that  these  are 
transmission  rates  at  the  sender  end.  Therefore,  the  transmission  rates  at  the  receiving 
node  are:  r,  = 20.6186*0.97  = 20,  = 20.202*0.99  = 20,  rj  = 15.4639*0.97  = 15 . 

6.5  Computational  Experience 

We  generated  243  small  problem  instances  with  given  ranges  in  Table  6-3.  We 
set  the  tightness  (the  ratio  of  total  task  size  to  total  bin  capacity)  to  50%  and  the  ratio  of 
number  of  time-fixed  tasks  to  number  of  size-fixed  tasks  to  1 . We  have  randomized  all 
other  factors  such  as  minimum  and  maximum  bandwidth  requirement  for  time-fixed 
tasks,  sizes  of  tasks,  bin  capacities  and  bin  costs  and  opportunity  costs  of  tasks. 


Table  6-3.  Factors  and  ranges  in  problem  setting 


Factor 

Range 

Number  of  bins  (|l|) 

3 to  8 

Number  of  tasks  per  bin 

2 to10 

Total  task  number  multiplier* 

1 to  HI-1 

* Total  number  of  tasks  = Total  task  number  multiplier  X Number  of  tasks  per  bin 

We  have  two  purposes  in  this  empirical  study.  First,  we  test  the  running  time  for 
P3  and  find  the  maximum  problem  size  that  can  be  solved  by  P3  in  a reasonable  time. 
For  this  reason,  the  number  of  QoS  levels  is  set  to  1 in  the  problem  setting.  Second 
purpose  is  to  compare  the  solution  quality  of  P2  and  P3  as  a lower  bound  for  PI . 
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GBD  is  a theoretically  elegant  procedure  and  compatible  with  P3.  However,  there 
are  computational  difficulties  related  to  the  problem  size.  GBD  works  well  for  small 
problems.  However,  when  the  problem  size  exceeds  a certain  level,  the  convergence 
rate  drops  drastically.  Table  6-4  shows  the  maximum  problem  sizes  where  GBD  works 
fast  for  us.  When  the  total  number  of  tasks  exceeds  15-20,  the  convergence  rate  starts 
slowing  down.  This  also  justifies  not  using  GBD  in  its  current  form  to  solve  PI  even 
though  in  section  6.4.2  we  showed  that  property  P technically  holds. 

Another  factor  that  affects  the  convergence  rate  is  the  number  of  tasks  per  bin. 
When  it  exceeds  3-4  convergence  rate  drops  drastically  as  given  in  Table  6-5.  The 
convergence  rate  drops  even  further  when  the  problem  size  increases.  For  example, 
when  there  are  10  bins  and  50  tasks  with  10  tasks  per  bin,  it  takes  more  than  8 hours  to 
complete  a single  iteration  after  the  first  several  iterations.  This  behavior  is  explained  by 
the  fact  that  step  (2)  of  GBD  for  this  problem  is  in  effect  a bin-packing  problem  with 
additional  constraints  on  the  binary  variables. 


Table  6-4.  Maximum  problem  sizes  where  GBD  works  fast 


Number  of  bins 

Total  number  of  tasks 

Number  of  tasks  per  bin 

3 

18 

9 

4 

16 

6 

5 

16 

4 

6 

15 

5 

7 

16 

4 

8 

15 

3 

Table  6-5.  Run  times  when  number  of  tasks  per  bin  is  10. 


Number  of  bins 

Total  number  of  tasks 

Run  time  (sec) 

6 

20 

0.812 

6 

30 

>900 

7 

20 

13.988 

7 

30 

>900 

8 

20 

3.034 

8 

30 

>1200 
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Looking  at  the  behavior  of  GBD  on  P3  we  observe  that  it  converges  to  a set  of 
resources  (bin)  quickly  and  spends  the  rest  of  the  time  swapping  tasks  among  these 
bins  while  recomputing  the  transmission  rates  accordingly  after  each  iteration.  Intuitively, 
finding  the  set  of  bins  should  take  a short  amount  of  time.  There  are  more  possibilities 
for  swapping  tasks  especially  when  the  transmission  rates  change  from  one  iteration  to 
the  other.  Depending  on  which  tasks  are  assigned  together  in  a bin  the  transmission  rate 
of  a fixed  task  might  change  from  one  iteration  to  another  just  because  the  associated 
dual  cost  of  the  bin  and  the  opportunity  costs  of  other  tasks  are  different.  When  the  total 
number  of  tasks  is  less  than  15-20,  we  expect  to  find  a solution  fast.  However,  it 
sometimes  takes  longer  to  find  a solution  since  there  may  be  more  swapping 
possibilities.  For  example,  when  there  are  8 bins  and  12  tasks  with  2 tasks  per  bin,  it 
takes  more  than  5 hours  and  101  iterations  to  terminate  the  program. 

Since  the  running  time  for  some  of  the  problems  in  the  dataset  is  too  large,  we 
selected  a subset  of  problems,  which  can  be  solved  in  reasonable  time.  We  further 
eliminated  some  of  the  problems  since  PI  could  not  be  solved  because  of  bandwidth 
violations  in  the  data  (even  though  P3  is  solvable  for  the  same  data  since  there  is 
enough  raw  capacity).This  left  120  problems  on  which  we  can  make  a comparative 
analysis  for  the  quality  of  two  lower  bound  procedures  for  PI . One  procedure  is 
described  in  this  chapter  (using  GBD  to  solve  P3)  and  the  other  in  Chapter  5 (Algorithm 
A1  used  to  solve  P2).  As  given  in  Table  6-6,  P3  yields  a tighter  lower  bound  for  PI  than 
P2.  We  were  able  to  show  that  Heuristic  A solves  36  of  the  120  problems  optimally 
(zpi=zp3).  In  10  instances,  GBD  was  stopped  prematurely  but  yielded  the  same  objective 
function  value  as  Heuristic  A. 


Table  6-6.  Quality  comparison  for  solutions  of  P2  and  P3 


Lower  Bound  Method 

P3 

P2 

Solution  Gap  between  Heuristic  A (%) 

7.65 

31.97 

Number  of  Problem  solved  Ootimallv 

36 

0 
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The  overall  average  gap  between  solutions  of  Heuristic  A and  P2  (solved  by 
procedure  A1)  is  reported  as  4.9%  in  Chapter  5.  Table  6-6  shows  that  this  gap  is  much 
larger  compared  to  the  results  of  Chapter  5.  We  believe  this  is  because  of  small 
problem  instances  and  low  tightness  levels  (at  50%).  Such  a low  tightness  inflates  the 
effects  of  prorating  (Vi  is  continuous  between  0 and  1 ). 

Thirty  six  out  of  120  problem  instances  were  solved  optimally  by  GBD  in  a 
reasonable  time  as  given  in  Table  6-7.  In  these  problem  instances,  the  gap  between  the 
solution  of  Heuristic  A (Za)  and  Lower  bound  procedure  A1  (Zai)  is  quite  large  because 
of  prorating. 

On  30  problems  GBD  procedure  did  not  terminate  after  15  minutes  and  were 
terminated  before  they  converged.  Their  solutions  are  given  in  Table  6-8.  Although  the 
GBD  procedure  is  terminated  before  convergence,  the  solution  quality  for  these  30 
problems  is  better  than  P2.  The  gap  between  LB  in  GBD  solutions  and  Heuristic  A 
solution  is  14.05%  while  the  gap  between  PI  and  P2  solution  is  20.66%.  This  is  because 
the  gap  between  UB  and  LB  during  GBD  is  closing  with  a decreasing  rate  and  the 
solution  close  to  optimal  could  be  obtained  in  early  iterations.  As  given  Figures  6-1  and 
6-2,  the  gap  decreases  very  fast  during  early  iterations.  Then  the  rate  is  decreasing 
during  later  iterations. 
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Figure  6-1.  LB  and  UB  values  during  GBD  for  problem  instance  “nB6_nJpB3_nTJ15" 
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Table  6-7.  Quality  comparison  for  problems  solved  optimally  in  reasonable  time. 


File  Name* 

Solution  of  GBD  (Zps) 

1-(Zp3/Zpi) 

1-(Zp2/Zpi) 

^P1 

LB 

UB 

Run  Time 

^P2 

nB3_nJpB5_nTJ5 

57.9196 

57.9196 

57.9196 

0.032 

21.075 

0 

63.61335 

nB3_nJpB7_nTJ7 

44.2702 

44.2702 

44.2702 

0.032 

19.3175 

0 

56.36455 

nB3_nJpB8_nTJ8 

47.7344 

47.7344 

47.7344 

0.063 

24.4971 

0 

48.68041 

nB3_nJpB9_nTJ9 

42.0191 

42.0191 

42.0191 

0.047 

20.4368 

0 

51.36307 

nB3_nJpB10_nTJ10 

46.8973 

46.8973 

46.8973 

0.094 

23.3156 

0 

50.28371 

nB4_nJpB2_nTJ2 

57.1882 

57.1882 

57.1882 

0.031 

4.47151 

0 

92.18106 

nB4_nJpB3_nTJ3 

42.7579 

42.7579 

42.7579 

0.041 

10.6246 

0 

75.15173 

nB4_nJpB4_nTJ4 

44.7764 

44.7764 

44.7764 

0.078 

18.5526 

0 

58.56612 

nB4_nJpB6_nTJ6 

40.7491 

40.7491 

40.7491 

0.109 

19.0519 

0 

53.24584 

nB4_nJpB6_nTJ12 

42.9685 

42.9685 

42.9685 

0.11 

40.3846 

0 

6.013475 

nB4_nJpB7_nTJ7 

46.7965 

46.7965 

46.7965 

0.063 

26.9691 

0 

42.36941 

nB4_nJpB9_nTJ9 

41.7341 

41.6153 

41.7341 

0.094 

21.5281 

0 

48.41604 

n B4_n J p B 1 0_nT J 1 0 

41.0404 

41.0404 

41.0404 

0.031 

22.6938 

0 

44.70376 

nB5_nJpB4_nTJ4 

42.3472 

42.3472 

42.3472 

0.062 

8.16515 

0 

80.71856 

nB5_nJpB5_nTJ5 

47.0073 

47.0073 

47.0073 

0.063 

22.4098 

0 

52.32698 

nB5_nJpB6_nTJ6 

58.0911 

58.0911 

58.0911 

0.156 

27.4357 

0 

52.77125 

nB5_nJpB7_nTJ7 

41.0336 

41.0336 

41.0336 

0.078 

17.7752 

0 

56.68135 

nB5_nJpB8_nTJ8 

48.5734 

48.5734 

48.5734 

0.281 

22.06 

0 

54.5842 

nB5_nJpB9_nTJ9 

49.7734 

49.7734 

49.7734 

0.125 

27.6991 

0 

44.34959 

nB5_nJpB10_nTJ10 

47.8555 

47.8555 

47.8555 

0.078 

27.8202 

0 

41.86624 

nB6_nJpB6_nTJ6 

45.409 

45.409 

45.409 

0.297 

30.2236 

0 

33.44139 

nB6_nJpB7_nTJ7 

56.4629 

56.4629 

56.4629 

0.078 

21.3587 

0 

62.17215 

nB6_nJpB7_nTJ14 

43.1532 

42.8502 

43.1532 

0.109 

31.3474 

0 

27.35788 

nB6_nJpB8_nTJ8 

40.2516 

40.2516 

40.2516 

0.157 

25.6894 

0 

36.17794 

nB6_nJpB8_nTJ24 

110.111 

94.1019 

110.111 

361.84 

83.2814 

0 

24.36596 

nB7_nJpB5_nTJ5 

41.4914 

41.4914 

41.4914 

0.125 

22.1095 

0 

46.71305 

nB7_nJpB5_nTJ10 

45.1488 

45.1488 

45.1488 

0.64 

39.9685 

0 

1 1 .47384 

nB7_nJpB8_nTJ8 

42.7518 

42.7518 

42.7518 

0.14 

19.5049 

0 

54.37642 

nB7_nJpB1 0_nTJ 1 0 

46.4042 

46.4042 

46.4042 

0.078 

24.149 

0 

47.95945 

nB8_nJpB3_nTJ3 

40.8293 

40.8293 

40.8293 

0.031 

6.51321 

0 

84.04771 

nB8_nJpB4_nTJ4 

41.1701 

41.1701 

41.1701 

0.234 

25.1045 

0 

39.02249 

nB8_nJpB5_nTJ5 

42.4382 

42.4382 

42.4382 

0.078 

1 1 .9273 

0 

71.8949 

nB8_nJpB6_nTJ6 

50.2975 

50.2975 

50.2975 

0.372 

29.1677 

0 

42.00964 

nB8_nJpB7_nTJ7 

41.8528 

41.8528 

41.8528 

0.109 

19.7688 

0 

52.76588 

nB8_nJpB7_nTJ14 

49.2502 

49.2502 

49.2502 

0.469 

45.8169 

0 

6.971139 

nB8  nJpB8  nTJ8 

43.2737 

43.2737 

43.2737 

0.141 

21.0976 

0 

51.24614 

Average 

0 

49.06241 

* In  file  name  nBx_nJpBy_nTJz,  x,  y and  z represents  number  of  bin,  number  of  jobs  per  bin  and 
number  of  total  jobs  respectively. 
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Figure  6-2.  LB  and  UB  values  during  GBD  for  problem  instance  “nB6_nJpB3_nTJ15" 


Table  6-8.  The  solution  quality  of  problem  instances  terminated  after  15  minutes. 


File  Name 

Zpi 

Solution  of  GBD  (Zps) 

LB  UB  Run  Time 

Zp2 

1-(Zp3/Zpi) 

1 -(Zp2/Zpi) 

nB3_nJpB9_nTJ18 

127.547 

122.074 

127.547 

1826.63 

72.749 

4.290967 

42.96299 

nB4_nJpB6_nTJ18 

109.409 

96.4359 

104.52 

1826.37 

70.3126 

11.85743 

35.73417 

nB4_nJpB7_nTJ21 

101.951 

97.7896 

101.951 

1894.33 

79.6049 

4.081765 

21.91847 

nB4_nJpB8_nTJ24 

85.1748 

81.4098 

85.1748 

1811.79 

62.3209 

4.420322 

26.83176 

nB4_nJpB9_nTJ27 

94.7015 

79.179 

92.3243 

2040 

87.1537 

16.39098 

7.970096 

nB4_nJpB10_nTJ20 

118.369 

106.985 

114.643 

1810.96 

74.9384 

9.617383 

36.69086 

nB4_nJpB10_nTJ30 

83.4437 

71.7877 

83.4437 

3684.7 

65.71.47 

13.9687 

21.24666 

nB5_nJpB4_nTJ16 

101.118 

95.7931 

101.118 

905.647 

81.5709 

5.266026 

19.33098 

nB5_nJpB5_nTJ20 

120.507 

112.948 

118.255 

901.647 

100.204 

6.272665 

16.84798 

nB5_nJpB6_nTJ18 

118.089 

114.447 

118.09 

1089.95 

70.7125 

3.084115 

40.11932 

nB5_nJpB6_nTJ24 

116.768 

90.1685 

114.854 

964.545 

104.192 

22.77979 

10.77007 

nB5_nJpB7_nTJ28 

140.029 

109.723 

121.07 

971.322 

122.346 

21.64266 

12.6281 

nB5_nJpB10_nTJ30 

108.774 

94.9014 

108.207 

961.421 

96.2886 

12.7536 

11.47829 

nB6_nJpB5_nTJ15 

122.734 

107.09 

115.992 

917.105 

87.0152 

12.74626 

29.10261 

nB6_nJpB5_nTJ20 

115.827 

107.063 

110.284 

911.426 

100.046 

7.566457 

13.62463 

nB6_nJpB6_nTJ18 

83.8347 

80.5136 

83.8347 

951.333 

57.0293 

3.961486 

31.97411 

nB6_nJpB7_nTJ21 

88.8609 

81.7144 

88.8609 

1243.91 

71.6969 

8.042345 

19.31558 

nB6_nJpB8_nTJ32 

107.011 

89.1413 

101.085 

966.95 

101.186 

16.69894 

5.443366 

nB7_nJpB4_nTJ12 

98.0832 

75.2191 

85.2256 

910.7 

67.8102 

23.31092 

30.86461 

nB7_nJpB4_nTJ16 

104.197 

90.7137 

99.8945 

903.587 

92.8148 

12.9402 

10.92373 

nB7_nJpB5_nTJ20 

113.944 

77.0325 

102.983 

962.428 

99.5346 

32.39442 

12.64604 

nB7_nJpB6_nTJ18 

99.5377 

82.1987 

96.5764 

1032.12 

73.3182 

17.41953 

26.34128 

nB7_nJpB7_nTJ21 

88.676 

85.3058 

88.676 

1742.65 

62.4304 

3.800577 

29.59719 

nB7_nJpB8_nTJ24 

100.518 

87.5492 

102.47 

924.803 

92.9894 

12.90197 

7.489803 

nB7_nJpB10_nTJ30 

87.8644 

73.8671 

85.2971 

920.727 

83.3081 

15.93057 

5.185604 

nB8_nJpB6_nTJ24 

113.227 

104.122 

111.611 

1319.01 

96.8851 

8.041368 

14.43286 

nB8_nJpB7_nTJ21 

112.33 

77.8772 

112.33 

982.822 

78.4212 

30.67106 

30.18677 

nB8_nJpB7_nTJ28 

95.4119 

66.1867 

94.2153 

1363.77 

85.0688 

30.63056 

10.84047 

nB8_nJpB9_nTJ27 

89.194 

72.815 

89.194 

995.646 

75.1499 

18.36334 

15.74557 

nB8  nJpBIO  nTJ30 

101.779 

71.5371 

92.1578 

1105.4 

79.8304 

29.7133 

21.56496 

Average 

14.05199 

20.6603 
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In  summary,  we  develop  a better  lower  bound  by  using  GBD.  We  show  that  its 
quality  is  much  better  than  that  of  P2  formulated  and  solved  in  Chapter  5.  We  show  that 
GBD  solves  small  problems  in  reasonable  time.  We  also  show  that  we  obtain  better 
solutions  compared  to  those  obtained  from  P2  even  if  we  stop  the  GBD  procedure  after 
reasonable  time  before  it  terminates. 


CHAPTER  7 

CAPACITY  ACQUISITION  AND  TASK  ALLOCATION  SUBJECT  TO  QOS 
REQUIREMENTS  WITH  TAX  - BAND  PRICING 

In  this  chapter,  we  formulate  a cost  minimization  problem  subject  to  QoS  and 
capacity  constraints.  We  consider  the  tax-band  pricing  scheme  suggested  by 
Courcoubetis  et  al.  (2000)  in  which  each  supplier  has  a convex  piecewise-linear  cost 
function  for  each  resource  offered  (Figure  7-1).  For  simplicity  we  assume  that  capacity  is 
unlimited  since  there  is  no  upper  limit  on  how  long  the  resource  is  available. 
Courcoubetis  et  al.  claim  that  tax-band  pricing  reduces  bursty  traffic  since  customers  are 
likely  to  reduce  such  demand  to  avoid  paying  more.  Suppliers  in  general  would  prefer 
having  more  customers  with  less  capacity  demand  rather  than  fewer  customers  with 
high  capacity  demand.  When  customers  demand  more  capacity  the  possibility  of 
network  congestion  increases  and  the  supplier  has  fewer  degrees  of  freedom  to 
maintain  QoS  levels.  Moreover,  a large  customer  pool  with  small  demand  should  also 
create  a more  uniform  network  usage.  Therefore,  suppliers  would  prefer  to  have  a large 
customer  pool  with  less  demand  per  customer  than  small  customer  pool  with  high 
demand  per  customer. 

Given  this  pricing  structure  the  customer  has  to  decide  how  much  capacity  to 
acquire  and  how  to  allocate  tasks.  As  in  Chapter  5,  for  tractability  we  relax  the  due  date 
constraints  and  assume  that  all  tasks  and  resources  are  available  at  time  zero  (Ei=0). 

We  also  assume  that  real-time  applications  have  desired  transmission  rates  and  any 
deviation  from  that  creates  an  opportunity  cost  for  the  customer.  In  spirit,  this  is  the 
same  objective  function  used  for  PI  that  trades  the  quality  cost  with  the  capacity  cost  but 
the  cost  structure  for  capacity  is  quite  different  in  this  problem. 
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In  this  chapter,  firstly,  we  present  tax-band  pricing.  Secondly,  a formulation  and  a 
solution  procedure  based  on  GBD  are  described.  Finally,  a more  realistic  formulation 
with  tax-band  pricing  considering  bandwidth  and  availability  of  resources  as  decision 
variables  instead  of  capacity  itself  is  given. 

7.1  Tax-Band  Charging  Scheme 

We  assume  that  the  supplier  i sells  capacity  (Si)  based  on  the  tax-band  pricing 
structure  (Figure  7-1).  In  tax-band  charging  scheme,  the  unit  cost  of  the  capacity  is  a 
step  function  where  slope  jumps  from  Cii  to  0,2  at  ^.because  of  diseconomies  of  scale. 
Such  a cost  function  should  discourage  demanding  large  capacity  from  a single  supplier. 


Figure  7-1.  Sample  price  curve  for  resource  i 

The  total  cost  is  defined  as  a function  of  acquired  capacity  as  given  below. 


c,|5,.  when  5,.  < 5,. 

<=n  (Si~Si)  + when  S.  > 


(7.1) 


Equation  (7.1)  is  a piecewise  linear  convex  function  and  we  assume  that  there 
are  two  different  unit  costs  (c,.,<  c,^ ). 

7.2  Problem  (P4)  Formulation  and  GBD  implementation  with  Tax-Band  Pricing 

/e/  ie/  j&Aj 
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subject  to 

Z ^kyik  + Z ^A^/jyy  ^ e / 

k^A^  jeAj 

(7.2) 

z,>c-Ai  yiel 

(7.3) 

+ yiel 

(7.4) 

V/G^,  Vfe/ 

(7.5) 

V/g/ 

^ Q^ij  vy  G j,  Vi  G / 

I yij  = L vy  € J 

i<=  J 

(7.6) 

yy  e {0,1} , z.,S^,Kj  >0  is  I,\fj  s J 

In  P4,  Si  is  the  acquired  capacity  of  resource  i and  z,  is  total  cost  of  that 
purchase.  Total  cost  is  constrained  by  equations  (7.3)  and  (7.4).  Either  (7.3)  or  (7.4)  is  a 
binding  constraint  depending  on  the  capacity  S|  selection. 

Let  u.,a.,b.,o)y,ii/y  islj  sA,.  be  the  dual  variables  corresponding  to 
constraint  sets  (7.2),  (7.3),  (7.4),  (7.5)  and  (7.6)  respectively.  Then,  the  RMP  of  GBD  for 
a given  fixed  y..  at  any  iteration  with  dual  solution  Ui,a.,^,ay,\py  can  be  written  as 


Max 

subject  to 


“Z  - Z Z + Z Z wlyu 

iel  iel  jsAj  iel 


i€l 


te/f,  jeAj 

^Z  Z ] - Z Z ] 

/e/  j^Aj  ieJ  JeAp 

+Z  ^ c,,  (5,  - 5;.)  - c,,5;. ) 


y^  <Max 


> 


(7.7) 
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Notice  that  some  of  the  terms  in  Inequality  (7.7)  are  independent  of  r. . We  can  then 
rewrite  (7.7)  as  follows. 


Z 

»€/ 


“1  \ 


- Z Z + Z - Z 'yijyi/j' + (^.-2  - c,-, ) 


JeAr 


keA, 


J^Aj- 


j^Aj- 


+Mca\ ^z.(a.  +b.-\)  + J]S,{u.a.  -a.c,,  -^,-2) 

' [/€/  le/ 

+ A/ox  X 0 Z [yij  (^1  - 

' -£/  J 


(7.8) 


Si  andz,.  depend  on  as  given  in  (7.2),  (7.3)  and  (7.4)  but  rj  can  be  set 


independent  of  any  other  variable.  Given  r.  and  y^,  S,  is  the  minimum  upper  bound  that 

makes  rj  feasible  and  Zj  is  the  associated  cost  of  doing  that.  If  the  term 

{uiUi  - afu  - ^Ci2 ) > 0 for  some  i in  (7.8)  then  5,.  ->  00  since  S|  has  no  upper  bound  and 

we  are  maximizing  in  (7.8).  This  however  does  not  make  sense  in  the  context  of  this 
problem  since  we  know  that  there  is  always  a finite  capacity  that  needs  to  be  purchased 

and  there  is  at  least  one  bounded  feasible  solution.  So  long  as  the  terms  (a).  + ^ - l) 

and  -a;.c,.,  -^c.^jare  both  negative  or  zero  it  is  possible  to  use  the  S|  and  Zi  values 
from  the  previous  solution  as  the  Ziand  Sj  values  that  maximizes 
Z z,  ^ - 1)  + Z - biCi2 ) since  they  are  the  lowest  values  possible  that 

le/  /'e/ 

are  feasible.  Theorem  6.1  is  still  valid  for  computing  the  ry  values. 

Note  that  S|  and  z,  values  that  are  substituted  from  the  previous  solution  is  just  an 
estimate.  It  is  possible  to  update  Si  and  Z|  of  iteration  t after  solving  SUB  (Step  3)  of 
iteration  t and  resubstituting  these  updated  values  in  the  t‘^  cut  while  creating  the  t+1*‘ 
cut.  Hence,  the  correct  form  of  a cut  is  created  in  the  next  iteration.  Using  this  procedure 
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caused  no  convergence  issues  in  our  computational  analysis  but  a rigorous  proof  of  the 
validity  of  this  cut  is  left  as  future  research. 

7.2.1  Special  Cases  for  P4 

7.2. 1.1  When  one  supplier  dominates 

If  c,.,  < c^,  yksi- {/} , yy  - 1,  Vj  (as  long  as  Qj  are  the  same)  since  resource  i 

is  the  cheapest  resource  for  all  range  of  capacity  if  the  QoS  constraints  allow  this 
assignment.  If  not,  the  problem  size  can  be  reduced  drastically  by  solving  P4  only  for 
those  tasks  that  could  be  assigned  to  resource  i. 

7.2.1.2  Task  splitting 

So  far  we  assumed  that  each  task  is  assigned  to  a single  resource.  Given  the 
right  setup  it  would  actually  be  possible  to  split  a task  and  submit  it  through  multiple 
networks.  Such  a formulation  (that  is  also  valid  in  the  context  of  PI)  simplifies  the 
problem  drastically.  Hence,  we  can  formulate  the  problem  as  follows. 


subject  to 


Z ^jyy  + Z ^i^jyu  ^ ^ ^ 


^ Qiyy  Vy  G J,  V/  G / 


S yy  = Vy  G J 

1C  I J 


Xj<Xj,  j&J 


z,.,5.>0,V/g7 
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where  = d^tj,  Xj  -Y/'Atj,  represents  the  fraction  of  job  j assigned  to  bin  i,  Xj  is 


the  amount  of  data  transmitted  and  c° 


It  is  easy  to  see  that  the  above  formulation  is  similar  to  P2  of  Chapter  5.  S-,  the 


amount  of  capacity  that  is  used  in  a given  resource,  is  equal  to  P^L.  in  P2.  z.,  the 
acquisition  cost  for  that  resource,  is  equal  to 

cf  — fj in  P2  since  ^ ^ ^i^jyy  - a, .5,.  and  ^ - c.,  where  k=1 

®/PiA-  ye-4,  j^Aj  PiA 

or  2. 

Algorithm  A1  of  Chapter  5 that  solves  P2  can  be  modified  and  used  for  solving 
P5.  While  assigning  tasks  whenever  we  split  a task  we  need  to  set  c.^  as  the  unit  cost  for 

a given  resource  and  resort  the  resources  in  ascending  order  of  their  unit  costs  and 
continue  assigning  tasks  with  the  resorted  resources.  When  crashing,  if  a resource  has 
been  allocated  tasks  whose  total  size  exceeds  5',. , crash  the  tasks  until  it  gets  smaller 

than  Sj.  At  this  point  set  c,.j  as  the  unit  cost  for  that  resource  and  resort  the  resources  in 

ascending  order  of  their  unit  costs  and  continue  reducing  size  of  tasks  with  resorted 
resources. 

To  get  an  idea  of  the  quality  of  the  solution  of  P4,  we  use  P5  as  a lower  bound 
since  it  is  a relaxation  of  P4. 

7.2.2  Computational  Experience 

For  our  empirical  study,  we  further  extend  the  problem  set  used  in  Chapter  6 by 
increasing  the  number  of  bins  from  8 to  10.  There  are  two  purposes  in  this  empirical 
study.  First,  we  test  the  running  time  for  P4  and  find  how  fast  GBD  is  working  with  the 
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tax-band  pricing  scheme.  Second  purpose  is  to  measure  the  solution  quality  of  P4  by 
comparing  the  result  with  P5. 

Although  the  capacity  of  each  resource  is  not  fixed,  P4  still  looks  like  a BPP.  The 
convergence  rate  for  the  problem  is  still  slow.  The  average  of  the  running  time  ratios  of 
P4  to  P3  considering  the  same  data  set  used  in  Chapter  6 is  less  than  0.38.  Table  7-1 
shows  that  there  is  a significant  running  time  improvement  in  some  problem  instances  of 
P4  when  compared  to  problem  instances  of  P3  of  the  same  size.  Even  though,  the 
running  time  for  P4  is  greater  than  P3  in  a few  problems,  the  difference  is  insignificant. 
Table  7-1 . Run  time  comparison  of  GBD  on  P3  and  P4 


File  Name 

Run  Time  (sec) 

Tp4  Tp3 

Tp4/  Tp3 

nB4  nJpB7  nTJ21.txt 

1.593 

1894.33 

0.000841 

nB5_nJpB3_nTJ12.txt 

4.583 

124.731 

0.03674 

nB5_nJpB6  nTJ24.txt 

2.010 

964.545 

0.002084 

nB5_nJpB8_nTJ32.txt 

214.936 

>900 

<0.238818 

nB6_nJpB3_nTJ12.txt 

18.979 

727.95 

0.026071 

nB6_nJpB6_nTJ30.txt 

14.937 

>1200 

<0.012448 

nB7_nJpB5_nTJ25.txt 

531.101 

>1200 

<0.4426 

nB7_nJpB8_nTJ40.txt 

33.354 

>1200 

<0.0278 

nB7_nJpB9_nTJ27.txt 

258.291 

>1200 

<0.2153 

NB8_nJpB5_nTJ25.txt 

614.43 

>1200 

<0.5120 

NB8  nJpB8  nTJ24.txt 

7.844 

>1500 

<0.005 

In  order  to  test  the  quality  of  solution  of  P4,  we  compared  its  results  with  the 
lower  bound  obtained  by  solving  P5.  We  tested  260  problem  instances  and  found  the 
average  ratio  of  the  solution  of  P4  to  the  lower  bound  as  1 .003012.  155  out  of  260 
problem  instances  are  solved  optimally.  Comparisons  of  some  of  the  problem  instances 
are  given  in  Table  7-2.  The  three  columns  under  P4  document  the  upper  bound,  lower 
bound  and  computational  time  of  GBD.  Column  Zps  reports  the  solution  of  P5  and  the 
next  two  columns  documents  the  difference  and  the  ratios. 

In  most  cases  the  solution  of  P4  is  either  equal  or  very  close  to  the  solution  of  P5 
since  the  problem  sizes  are  rather  small.  The  difference  occurs  due  to  a single  task  per 
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resource  where  the  solution  procedure  has  to  decide  if  the  last  job  assignment  uses 
more  capacity  than  The  formulation  of  P5  allows  for  splitting  the  task  whereas  P4 


does  not.  It  either  incurs  the  extra  cost  due  to  c.^  or  assigns  it  to  another  resource 

creating  the  same  situation  for  another  task.  If  the  problem  size  increases,  there  could 
be  more  splitting  possibility  and  therefore  the  gap  between  the  solution  of  P4  and  P5 
could  be  larger. 

Table  7-2.  The  solution  comparison  of  GBD  and  LB  procedure  for  tax-band  scheme 


File  Name 

LB 

P4 

UB 

Run  Time 

Zp5 

Zp4  — Zp5 

Zp4  / Zp5 

nB4_nJpB4_nTJ8 

61.3146 

61.8165 

2.297 

61.7638 

0.0527 

1 .000853 

nB4_nJpB8_nTJ16 

48.5461 

48.9631 

2.766 

48.8199 

0.1432 

1 .002933 

nB5_nJpB3_nTJ12 

81.4477 

81.7233 

6.874 

79.6893 

2.034 

1 .025524 

nB5_nJpB8_nTJ16 

62.9558 

63.577 

42.281 

63.5714 

0.0056 

1.000088 

nB6_nJpB2_nTJ10 

115.321 

116.317 

41.512 

115.898 

0.419 

1.003615 

nB6_nJpB2_nTJ6 

99.0299 

99.7118 

0.86 

95.8681 

3.8437 

1 .040094 

nB6_nJpB3_nTJ12 

89.6828 

90.5813 

28.468 

89.5774 

1 .0039 

1.011207 

nB6_nJpB3_nTJ15 

110.617 

1 1 1 .694 

20.389 

109.237 

2.457 

1 .022492 

nB7_nJpB3_nTJ9 

93.9727 

94.4222 

6.578 

90.0573 

4.3649 

1 .048468 

nB7_nJpB7_nTJ21 

61.8673 

62.4362 

4.172 

62.4304 

0.0058 

1 .000093 

nB8_nJpB2_nTJ12 

146.379 

147.804 

11.718 

147.62 

0.184 

1.001246 

nB8_nJpB2_nTJ8 

80.6819 

80.7936 

2.812 

80.6514 

0.1422 

1.001763 

nB8_nJpB3_nTJ18 

92.7146 

93.5879 

28.046 

93.421 

0.1669 

1.001787 

nB8_nJpB4_nTJ16 

95.3594 

96.3135 

95.906 

96.3132 

0.0003 

1 .000003 

nB9_nJpB10_nTJ20 

57.6175 

58.1824 

56.786 

58.1823 

0.0001 

1 .000002 

nB10_nJpB2_nTJ8 

97.0671 

98.0277 

10.968 

97.7347 

0.293 

1 .002998 

nB10_nJpB7_nTJ14 

52.8047 

53.2988 

2.125 

53.2855 

0.0133 

1 .00025 

nB4_nJpB9_nTJ27 

79.8857 

86.7727 

920.191 

86.4736 

0.2991 

1 .003459 

nB5_nJpB10_nTJ40 

110.61 

117.615 

914.013 

117.397 

0.218 

1.001857 

nB5_nJpB5_nTJ20 

98.0389 

100.719 

911.06 

99.5541 

1.1649 

1.011701 

nB5_nJpB7_nTJ28 

113.007 

121.533 

947.659 

121.225 

0.308 

1.002541 

nB5_nJpB9_nTJ36 

106.754 

110.829 

955.056 

109.726 

1.103 

1.010052 

nB6  nJpB4  nTJ20 

123.889 

132.036 

910.448 

124 

8.036 

1 .064806 

Average 

4.3649 

1.009613 
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7.3  Summat7  and  Conclusions 

GBD  works  much  faster  in  with  the  tax-band  scheme  compared  with  all-you-can- 
send  pricing.  Even  stopping  it  prematurely  yields  very  good  results.  For  future  research, 
a more  realistic  tax-band  formulation  is  given  in  the  following  section. 

7.3.1  A More  Realistic  Formulation 

A more  realistic  formulation  for  tax-band  pricing  can  be  obtained  by  considering 
the  bandwidth  and  duration  of  resources  as  decision  variables.  In  this  case,  the  total  cost 
of  acquisition  is  a piecewise  linear  convex  function  of  bandwidth  as  given  Figure  7-2. 

The  unit  cost  of  bandwidth  after  p,.  increases  (C;"'  > cf' ) because  of  diseconomies  of 

scale. 

Total  Cost 


Figure  7-2.  Tax-band  charging  scheme  for  a given  QoS  level  for  the  supplier 
The  acquisition  cost  for  a given  resource  is  calculated  as  follows. 

Acquisition  Cost=  c]L.+  f.{p.)  (7.9) 


m)= 


if 


(7.10) 


The  acquisition  cost  has  two  terms.  The  first  term  is  related  to  the  duration  of  the 
resource  and  is  called  the  connection  cost.  The  second  term,  (7.10),  is  a function  of  the 
bandwidth  purchased.  There  is  a trade-off  between  these  two  cost  terms.  When  a 
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resource  with  high  bandwidth  is  purchased,  the  duration  required  to  complete  a fixed 
amount  of  jobs  decreases,  therefore  the  first  term  decreases  while  the  second  term 
increases.  Conversely,  purchasing  a low  bandwidth  resource  reduces  the  bandwidth 
related  costs  but  increases  connection  costs.  The  problem  formulation  is  given  in  the 
following  section. 

7. 3. 1.1  Problem  Reformulation 

(P6)  min  Y,  V,-  (c'A  + ^/ ) + Z Z (7- 1 ■' ) 

iel  iel  j^A-p 

subject  to 

Z + Z ^ Vf  € / (7.12) 

je/),  j&Aj 

z-  > c['  p.  \fi  el 

+ Vf€/ 

V/'G^,  V/e/ 

^Ai-y/j'  y/'e4>  Vie/ 

9 >0,:v,  VjeJ,  Vi  el 

Zt(,=i  y/e*/ 

iel 

yij<Vi  VieI,VjeJ 
rj,z.,/3.,L.>0  VieI,VjeA, 

P6  is  not  a trivial  problem  to  solve.  The  main  difference  between  P6  and  P4  is 
that  the  variable  in  P4  is  replaced  by  in  equation  (7.12)  and  the  acquisition  cost 
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in  (7.11)  depends  both  on  and  L|.  Therefore,  if  we  want  to  use  GBD  for  this  problem 
we  need  to  solve  a convex  minimization  problem  (rather  than  an  LP)  for  fixed  v.,yy . 

Figure  7-3  shows  that  the  solution  to  the  problem  of  finding  P;  and  L.  can  be 
found  by  comparing  the  slope  of  the  cost  function  to  the  derivative  of  the  capacity 
function  given  by  5 = P^*L. . The  slope  of  the  curve  at  p,.  is  given  below. 

5A-  " PI 


A Li 


PiP=S 


P,  P,- 

Figure  7-3.  Cost  and  capacity  curves  for  a given  capacity  for  the  supplier 

3z 

We  can  determine  by  checking  the  slope  of  the  cost  function  . If  it  is  greater 


than  , it  implies  /3.  < (5. . Extending  this  idea  we  can  find  p.,L.  for  a given  capacity  S 

Pi 


as  follows. 
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If 


c.'  ^ c] 

— <-=T<  — 

< Pf  c\ 


If  =>A=A  A 


A 


If 


Even  if  we  can  find  P^,L.  for  a given  capacity,  P6  is  not  trivial  to  solve  since  finding  the 
right  amount  of  capacity  for  each  resource  cannot  be  determined  easily  for  a given 
v.,yy , since  it  also  depends  on  Vj  as  given  in  (7.12).  We  leave  the  solution  procedure 
for  P6  as  future  research. 


CHAPTER  8 

CONCLUSION  AND  FUTURE  RESEARCH 

In  this  section  we  first  present  our  conclusions  and  summary  followed  by  a 
discussion  of  some  future  research  ideas. 

8.1  Conclusion  and  Summary 

This  study  investigates  a vendor  selection  and  resource  allocation  problem  in 
telecommunication  networks.  As  a rational  decision  maker  the  firm  has  to  determine 
which  network  resources  should  be  acquired  and  how  existing  and  scheduled  future 
tasks  should  be  assigned  to  these  resources.  In  the  process  the  customer  has  to  solve  a 
cost  minimization  problem  that  models  a trade  off  between  the  cost  of  acquiring 
resources  and  the  opportunity  cost  of  degradation  in  realized  quality  of  the  real  time 
applications.  The  formulations  consider  QoS  and  capacity  constraints  and  two  different 
pricing  schemes.  The  all-you-can-send  pricing  assumes  that  the  suppliers  charge  a 
fixed  price  for  a specific  bandwidth  and  duration,  whereas  the  tax-band  pricing  considers 
a piecewise  linear  convex  cost  function.  We  demonstrate  that  the  task  set,  which  is 
heterogeneous  in  nature,  presents  some  interesting  challenges.  We  also  suggest  explicit 
mechanisms  to  handle  the  QoS  requirements,  such  as  the  transmission  efficiency. 

We  show  that  the  resulting  problem  with  all-you-can  send  pricing  (PI)  is  at  least 
as  hard  as  2D-BPP  and  offer  a simple  heuristic,  that  can  take  advantage  of  the  trade-off 
between  quality  related  costs  and  capacity  costs,  since  obtaining  an  exact  solution  is 
relatively  hard.  We  determined  how  the  customer  behaves  under  different  pricing 
settings.  As  expected,  a customer  first  tries  to  select  the  cheapest  resources.  However, 
since  the  problem  is  a constrained  optimization  problem,  the  optimal  solution  requires 
switching  suppliers  because  of  delay  and  jitter  requirements  when  a single  supplier 
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could  not  satisfy  the  customer  concerning  QoS  even  if  it  is  the  cheapest  supplier. 
Moreover,  the  combinatoric  nature  of  the  problem  might  sometimes  force  the  customer 
to  switch  a supplier  and  select  more  expensive  resources  since  using  a more  expensive 
resource  with  a larger  capacity  might  reduce  wasted  capacity  decreasing  the  per  unit 
cost.  Most  models  discussed  in  Chapter  3 may  not  be  able  to  capture  this  behavior 
since  they  ignore  some  of  the  QoS  and  capacity  constraints. 

We  then  develop  two  lower  bound  formulations.  First  is  an  LP  relaxation  (P2)  and 
the  second  is  a relaxation  of  the  bandwidth  constraint  only  (P3).  The  LP  relaxation  is 
relatively  tight  when  resource  prices  are  set  randomly.  P3  is  a tighter  relaxation  and  is 
solvable  by  GBD,  an  exact  algorithm,  developed  by  (Geoffrion,  1972).  GBD  runs  fast  on 
small  problems  but  it  slows  considerably  on  large  problems  since  it  needs  to  solve  a 
BPP  at  each  iteration.  In  our  experience  when  the  problem  size  exceeds  15-20  tasks  or 
tasks  per  bin  count  exceeds  3-4,  the  convergence  rate  drops  drastically.  The  average 
gap  between  the  lower  bound  of  GBD  solutions  and  Heuristic  A solutions  is  14.05% 
when  it  is  prematurely  terminated.  Moreover,  we  may  obtain  a better  feasible  solution  to 
PI  by  using  GBD  rather  than  Heuristic  A by  stopping  the  GBD  procedure  after  a 
reasonable  time  since  we  have  observed  that  the  solution  quality  of  GBD  in  this  case  is 
still  quite  good. 

Finally,  we  show  that  P4  is  solved  much  faster  than  P3.  The  problem  structure 
changes  in  such  a way  that  the  packing  constraints  turn  into  soft  constraints  (variable 
right-hand-side  vs.  data).  The  results  show  that  if  one  supplier  offers  the  cheapest 
resource  in  all  range  of  capacity  in  P4,  the  customer  does  not  switch  (with  the  exception 
of  QoS  constraints).  In  general,  if  there  is  no  setup  cost  and  if  every  supplier  offers  tax- 
band  pricing,  an  optimal  strategy  for  the  customer  would  be  to  use  more  suppliers  so 
that  the  resource  usage  is  always  below  the  threshold  S.. 
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It  is  interesting  to  note  that  a more  practical  formulation  that  could  be  solved 
optimally  very  fast  is  the  case  where  we  can  assume  that  a given  task  can  be  split  into 
pieces  and  transmitted  over  multiple  resources  if  needed  (P2).  The  technology  to  do  this 
exists  today. 

8.2  Future  Research  Ideas 

Future  research  will  proceed  in  various  dimensions.  First,  we  will  investigate  in 
detail  whether  GBD  creates  better  feasible  solutions  than  Heuristic  A for  PI  when 
stopped  prematurely.  Moreover,  when  task  splitting  is  allowed,  PI  can  be  solved  in  a 
reasonable  time.  An  adaptation  of  Algorithm  A1  of  P2  to  two  dimensions  can  be  used  for 
this  purpose. 

Secondly,  we  will  focus  on  P6.  It  is  a very  realistic  formulation  of  the  original 
problem  and  it  presents  unique  theoretical  challenges. 

A possible  approach  to  get  around  the  “hard”  packing  constraints  of  PI  and  P3  is 
to  use  the  tax-band  formulation  as  a relaxation  of  the  all-you-can-send  pricing 
formulation.  When  c,.^  = oo , P4  resembles  P3.  Since  P4  is  relatively  easier  to  solve,  we 

might  be  able  use  a similar  procedure  to  solve  P3  by  first  converting  it  to  P4. 

A third  venue  for  future  research  is  incorporating  QoS  using  the  amount  of  delay 
and  jitter  rather  than  considering  them  as  simple  levels.  A probabilistic  quality  measure 
may  be  introduced  for  the  delay  and  jitter  components  of  the  QoS  requirement. 

Assuming  that  delay  has  a distribution  F we  can  form  constraints  such  as  \-F{d)  > p 

which  models  that  the  supplier  will  guarantee  a minimum  required  delay  d with 
probability  p . Finally,  supplier  switching  costs  and  penalty  costs  for  not  satisfying  QoS 
requirements  and  mixture  of  charging  schemes  can  be  added  to  the  model. 

It  would  be  relatively  easy  to  extend  this  model  for  a network  intermediary  by 
reformulating  it  as  a profit  maximization  problem  and  modifying  Constraint  4.11. 
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